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Final Report 


This report presents a synopsis of the research work performed under grant number NAG3-1507. 

The goals of the research effort, as described in the grant proposal, have been fulfilled. Specific 

accomplishments are itemized below. 

1. Experimental facilities have been developed. This includes a magnetic bearing test rig 
and an auxiliary bearing test rig. In addition, components have been designed, 
constructed, and tested for use with a rotordynamics test rig located at NASA Lewis 
Research Center. Photographs of the specific articles are shown in Appendix A. 

2. The tasks performed under Grant No. NAG3-157 are detailed below. Papers that 
document the accomplishments of the research effort are included in Appendix B. 

i. A study of the rotordynamics of an auxiliary bearing supported T-501 engine 
model was performed. A paper documenting this work was presented at the 
Symposium on Nonlinear and Stochastic Dynamics , held at the 1994 ASME Winter 
Annual Meeting, November 13-18, 1994, Chicago, Illinois. The paper has also 
been submitted to the ASME Journal of Vibration and Acoustics. An additional 
paper that discusses some of the work using the T-501 engine model has been 
presented at the First Industry /University Symposium on High Speed Civil 
Transport Vehicles, December 4-6, 1994, Greensboro, North Carolina. 

ii. An experimental/simulation study of auxiliary bearing rotordynamics has been 
performed. A paper that documents this work was presented at the 40 h ASME 
International Gas Turbine and Aeroengine Conference, June 5-8, 1995, Houston, 
Texas. An additional paper will be presented at the ASME International Mechanical 
Engineering Conference and Exposition, Nov. 17-22, 1996, Atlanta, Georgia. 

iii. A rotordynamical model for a magnetic bearing supported rotor system, including 
auxiliary bearing effects has been developed and simulation studies performed. A 
paper that documents this work was presented at the ASME 15th Biennial 
Conference on Vibration and Sound, Boston, Massachusetts, Sept. 17-21, 1995. 

iv. A finite element model for a foil bearing has been developed and studies of a rotor 
supported by foil bearings have been performed. A paper that documents the work 
was presented at the 1996 SDM Conference, held in Salt Lake City, Utah, April 
15-17, 1996. 

3. Two students affiliated with this project have graduated with M.S. degrees - April M. Free 
and James L. Lawen, Jr. Ms Free is currently employed by Pratt and Whitney Aircraft 
Engines, Inc. (West Palm Beach, Florida) and Mr. Lawen is employed by the General 
Electric Corporate Research and Development Center (Schenectady, New York). 
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Figure 1: Magnetic Bearing Supported Rotor Test Ri 
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Figure 2: Auxiliary Bearing Rotor Test Rig 
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Figure 3.b: Auxiliary Bearing Test Component for use with 

NASA/Lewis Rotordynamics Test Rig 
(Close-up View) 
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G. T. Flowers 
Fangsheng Wu 


Disk/Shaft Vibration Induced^o^ 
by Bearing Clearance Effects: 
Analysis and Experiment 


Department of Mechanical Engineering, 
Auburn University, 
Auburn, AL 


This study presents an investigation of the dynamics of a rotor system with bearing 
clearance. Of particular interest is the influence of such effects on coupled disk/ 
shaft vibration. Experimental results for a rotor system with a flexible disk are 
presented and compared to predictions from a simulation model. Some insights and 
conclusions are obtained with regard to the conditions under which such vibration 
may be significant. 


Background and Motivation 

Rotor systems typically consist of a shaft with one or more 
bladed disks attached. Disk and blade vibration are issues of 
important concern for stress analysts. Complex finite element 
models are developed to assess the dynamic stresses in such 
components and insure the integrity of the design. However, 
models for rotordynamical analyses are typically developed as- 
suming that disk flexibility effects are negligible. It is generally 
presumed that the disk is sufficiently rigid so as to not signifi- 
cantly impact rotor vibration over the operating region. 

There is a fairly large body of work documented in the litera- 
ture concerning studies of disk flexibility on rotors and turboma- 
chinery. A discussion of this previous work with regard to linear 
rotor systems is presented in Flowers and Ryan (1991). An 
excellent source for a comprehensive review of work in this 
area is presented by Davis ( 1989). 

Previous investigators into the area of coupled rotor/ disk 
vibration have noted that disk flexibility has little effect on 
critical speeds but that it may significantly influence higher 
natural frequencies of the rotor system (Chivens and Nelson, 
1975). One can draw the conclusion from these results that 
synchronous vibration due to imbalance will be little affected by 
disk flexibility. However, there are sources of higher frequency 
excitation that could serve to excite these natural frequencies. 
Perhaps the most obvious are multi- synchronous effects corre- 
sponding to the blade pass frequency (from fluid forces imping- 
ing on the rotor blades). Another potential source is from non- 
linearities that may be present in the system. For example, 
nonlinear bearing forces due to clearance effects may result in 
supersynchronous rotor vibration. There are quite a number of 
studies in the literature concerning the effects of bearing clear- 
ance (and the related phenomena of rubbing) on rotordynamical 
behavior. Some of the works that have most influenced the 
current study were conducted by Johnson (1962), Black 
(1968), Ehrich (1966), and Childs (1978), and Muszynska 
(1984). 

Advanced designs for many types of high speed rotating 
machinery that use magnetic bearings for support have been 
proposed and are in the development and construction stages. 
There are a number of such systems already in commercial 
use. Rotor systems supported by magnetic bearings must have 
backup bearings to provide support under overload conditions 
or if the magnetic bearing fails. Backup bearings are character- 
ized by a clearance between the rotor and the bearing such that 
contact does not occur under normal operating conditions. As 
a result, issues related to the effects of bearing clearance on 
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rotordynamical behavior are of current concern. The objective 
of the present work is to develop an understanding of possible 
coupling between the dynamics of disk and shaft that may be 
induced by bearing clearance effects and to provide guidance 
to designers concerned with such systems. 

Experimental Model 

In order to investigate whether bearing clearance can lead to 
significant coupling between rotor and disk vibration, experi- 
mental tests were performed with a rotor test rig. A drawing of 
the test rig is shown in Fig. 1. 

The rotor used in this study has two basic components: a 
flexible disk and a shaft. The steel shaft is 0.375 in. diameter 
and 16.0 inches long. The flexible disk is a circular aluminum 
plate 0.0125 inches thick and 14.0 inches in diameter. The 
natural frequency of its lowest one nodal diameter bending 
mode is about 40 Hz. It is attached with four bolts to a 2.0 inch 
diameter hub mounted with a press fit to the center of the shaft. 

The rotor shaft is supported by bushings at two ends. The 
left end of the shaft is placed directly in the bushing base. This 
support provides a force but only a minimal couple so that the 
shaft is effectively free to pivot about this point. The right 
bearing is supported by a special device, which is designed to 
provide stiffness and lo simulate bearing clearance. A diagram 
for this device is shown in Fig. 2. It has two sets of springs. 
The softer springs are to support the weight of the rotor, and 
the harder springs act as the main stiffness for the system. The 
stiffness in the x and y directions of each spring set may be 
different. The clearance is adjustable by turning the inner ta- 
pered ring in or out. In this experiment, the spring constant of 
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the soft spring is 1.9 lb/ in. The spring constant of the hard 
spring is 8.1 lb/in. in the x, (horizontal) direction and 19.7 lb/ 
in. in the x 2 (vertical) direction. 

The rotor is driven by a speed adjustable motor. The speed 
of the rotor can be controlled by turning the speed control 
knob on the front panel of the control box. Shaft vibration was 
measured using eddy current proximity displacement sensors 
fixed so as to measure displacement in the vertical and hori- 
zontal directions. Strain gages mounted to the disk were used 
to measure the disk vibration. The strain gage signals were 
transmitted from the rotating disk to the signal conditioner using 
a set of slip rings affixed to the shaft. The displacement signals 
of the shaft and disk were recorded and analyzed with a signal 
analyzer. 

Simulation Model 

The simulation model is similar to that developed in an earlier 
study (Flowers and Wu, 1993). The primary differences are 
that rotational stiffening of the disk have been taken into account 
and that damping is included in the bearing support forces. A 
schematic diagram of the model is shown in Fig. 3. It consists 
of a rigid shaft, a rigid hub, a flexible disk and a support with 
a symmetric clearance. The equation development is based on 
the following considerations: 

( 1 ) The disk is assumed to flex only in the lateral direction. 

(2) Only rotational vibration is considered. This is because 
only the rotational motion is coupled with the lateral disk 
vibration. One nodal diameter disk vibration is assumed. 



( 3 ) The rotor speed is constant. 

After some mathematical manipulations, the dimensionless 
equations of motion for the whole system can be obtained as 

x, + r,x 2 + r 2 x 3 + 2r 2 x 4 + /,, = cos r (1) 

x 2 - r,x, + r 2 x 4 — 2r z x% + f n2 = sin r (2) 

*3 + 2 x 4 + l&UJiXi + (lu'i - l)x 3 = -(if, + 2*2) ( 3 ) 

x 4 - 2X7, + 2£ 3 u/ 3 x 4 + - 1 )X A = -(*2 - 2x y ) (4) 


fni = ~<t> 

[0 + 2 * ,w,x '] 

fm = ~4> 

[{ l 7^T7l) ( " 2X2) + 2 ^ 2Xl ] 


6 = * 


u — 

e 


— 1 if Vx 1 + x 2 > S 


0 otherwise 


In the above equations, x,, x 2 , x 3 and x 4 are scaled by e, the 
imbalance eccentricity. Time is scaled by the rotor speed, H. 
f „ , and/ n2 are the nonlinear moments due to contacting between 
the rotor and the bearings. 

The effective stiffness of a spinning disk can be strongly 
influenced by the spin speed. Based upon earlier work ( Wu and 
Flowers, 1992), the lowest natural frequency of such a disk 
can be written as 


2 19 + v u ! 30 

where o; 30 is the natural frequency of the nonspinning disk. This 
expression is used in Eqs. (3) and (4) to account for spin 
stiffening effects. 


Nomenclature 


e = imbalance eccentricity 
fn 1 » fni = nonlinear bearing forces 

r, =s gyroscopic mass influence ratio 
r 2 = disk mass influence ratio 
x, , x 2 = shaft degrees of freedom 
x 3 , x 4 = disk degrees of freedom 


A = bearing clearance 
f) = rotor speed 

a;, = natural frequency of rotor support 
in x, direction 

c o 2 — natural frequency of rotor support 
in x 2 direction 

cj 3 = natural frequency of rotor disk 


r = fit 

= damping ratio of rotor support in 
x, direction 

f 2 = damping ratio of rotor support in 
x 2 direction 

£3 = damping ratio of rotor disk 

C) = ± 

dr 
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Comparison of Experiment and Simulation 

A number of tests were conducted using the experimental 
setup described in the previous section. In addition, studies for 
differing parametric configurations have been conducted using 
the simulation model described above. The basic results are 
very similar to those obtained by Flowers and Wu ( 1993). The 
results presented here are typical. The following discussion is 
directed at comparing the predictions of the simulation model 
with experimentally observed responses, with the objective of 
obtaining insight into the behavior of flexible disk rotor systems 
with bearing clearance effects. 

First, the linear characteristics of the system were examined. 
Figure 4 shows experimental data of the natural frequencies for 
the coupled disk/ shaft vibration mode as a function of rotor 
speed obtained using the impulse excitation method. Using the 
data from such tests and additional measurements and calcula- 
tions, the linear stiffness, damping, and mass characteristics of 
the experimental rotor system were identified. The numerical 
values are shown in Table 1. 

Figure 4 also shows the natural frequencies as a function of 
rotor speed obtained from the simulation model. These results 
agree relatively well with the experimentally obtained data and 
serve to validate the structural parameters selected for the simu- 
lation model. 0, 3fi, and 5Q lines are also shown on the figure. 
The intersections of these lines with the natural frequency curve 
provides significant insight into the nonlinear behavior of the 
flexible disk rotor system, as is discussed below. 

Next, the nonlinear behavior of the system was investigated. 
A clearance value was selected and the rotor speed was slowly 
increased over a range of rotor speeds, with the resulting re- 
sponse amplitudes and frequencies measured and recorded. 

It is important to note that many vibration frequencies are 
possible. During the course of the study, lf2 — 60 frequency 
components were observed for certain parametric configura- 
tions. In addition, frequency components that are fractions of 
the rotor speed were observed for certain speed ranges and 
initial conditions. However, for a rotor with a symmetric clear- 
ance and minimal frictional effects, the odd integer multiples 
of the rotor speed appear to generally be the primary supersyn- 
chronous frequency components. The current discussion will 
concentrate on those frequencies. However, the basic conclu- 


sions should be directly applicable to other supersynchronous 
vibration frequencies as well. 

Figure 5 shows the experimentally measured amplitudes of 
the Q. 3H, and 5f2 components of the shaft response for various 
rotor speeds. Simulation results are also shown on Fig. 5. These 
results were obtained using the harmonic balance method and 
were verified at selected points through direct numerical integra- 
tion of the equations of motion. There is relatively good 
agreement between the simulation results and the corresponding 
experimental data. Certainly the basic trends match quite well. 

An important observation that can be made from examination 
of these figures is that the peaks of the respective components 
correspond to the intersections of the Q, 3f L and 5fl lines with 
the natural frequency curves. It appears that the nonlinear effects 
are serving to excite the coupled modes of the rotor/disk system. 
This is true for both the forward and backward whirl modes of 
the rotor system. 

The peaks occurring at a rotor speed of 0.5 1 correspond to 
the intersections of the 30 line with the natural frequency curve. 
They relate to a backward whirl mode. The peaks occurring at 
rotor speeds of 0.33 and 0.485 correspond to the intersections 
of the 50 line with the natural frequency curve. Those occurring 
at 0.33 relate to a backward whirl mode and those occurring at 
0.485 relate to a forward whirl mode. Note that the amplitudes 
of the supersynchronous components associated with the for- 
ward whirl mode of the disk/shaft vibration are much lower 
than those associated with the backward whirl mode. This result 
was observed for ail the parametric configurations examined in 
the course of this study. 

The Q, 3H, and 5 Q components contribute (more or less 
significantly ) to the overall response at the corresponding inter- 
section point. However, it is important to note that the Q compo- 
nent is increased rather dramatically at the intersection points 
as a result of the presence of disk flexibility. Apparently, the 
nonlinear coupling between the various frequency components 
serves to produce this effect. 

Conclusions 

A study of disk/ shaft vibration induced by bearing clearance 
effects has been presented. Both experimental tests and simula- 
tion results have been presented. The responses predicted by 
the simulation model and those observed experimentally agree 
quite well. The behavior of both the experimental test rig and the 
simulation model was quite sensitive to changes in parametric 
configuration. As a result it is quite difficult to predict exactly 
how a certain system is going to behave in the presence of a 
bearing clearance effect, i.e., whether or not coupled disk/ shaft 
vibration will be significantly excited. However, a few general 
conclusions and guidelines can be drawn from this study. 

1. An understanding of the mechanism for coupled disk/ 
shaft vibration induced by clearance effects has been 
obtained. The nonlinear effects have served to produce 


Table 1 Simulation model parameters 


Parameter 

Value 

ri 

0.5 

r 2 

0.1 

6 

0.06 

6 

0.06 

6 

0.01 

Wi 

0.55 

U>2 

0.85 

U>3 

1.8 

V 

0.3 

5 

1.2 
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Fig. 5(a) Q component of shaft vibration response Fig. 5(b) 3fi component of shaft vibration response 

Simulation Model 



superharmonics that excite coupled disk shaft modes. 
Conversely, the additional degree of freedom provided 
by the disk flexibility has served to exaggerate the fre- 
quency components of the response, resulting in higher 
amplitudes at the rotor speeds corresponding to the re- 
spective natural frequency /multi-synchronous line inter- 
section points. 

2. Bearing clearance effects primarily serve to excite back- 
ward whirl modes for coupled disk/ shaft vibration. For- 
ward whirl modes may be excited but the associated 
amplitudes appear to be much lower. 

3. The bearing support stiffnesses used in this study were 
not symmetric. These nonsymmetries, together with the 
clearance, have served to produce the supersynchronous 
vibration that excites the coupled disk/shaft vibration 
modes. Please note that the influence of nonsymmetrical 


stiffnesses and bearing clearance in this regard has been 
discussed in further detail by Flowers and Wu ( 1993). 

4. From a design perspective, the Campbell diagram is a 
useful tool to predict when such behavior may occur. 
The intersections of the corresponding supersynchronous 
line with the coupled disk/shaft vibration mode fre- 
quency curve will indicate at what rotor speeds peak 
responses for that frequency component are likely to oc- 
cur. Whether or not such behavior will actually occur 
and the relative significance of such effects depends very 
strongly on the imbalance, clearance, damping, and stiff- 
ness values. 

5. It appears to be relatively difficult to excite coupled disk/ 
shaft vibration with bearing clearance effects. The occur- 
rence of such behavior was very sensitive to rotor speed, 
as judged by the relative difficulty encountered during 


Journal of Vibration and Acoustics 


APRIL 1996, Vol. 118/20? 



the experimental work in obtaining the peak amplitude 
responses. This conclusion is verified by the simulation 
work that indicates that the behavior occurs only over a 
very limited range of rotor speed. The supersynchronous 
frequency must almost exactly coincide with a natural 
frequency in order to excite the behavior. 
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DYNAMIC BEHAVIOR OF A MAGNETIC BEARING 
SUPPORTED JET ENGINE ROTOR WITH AUXILIARY BEARINGS 
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H. Xie 
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Auburn University 
Auburn, Alabama 


ABSTRACT 

This paper presents a study of the dynamic behavior of a rotor system supported by 
auxiliary bearings. The steady-state behavior of a simulation model based upon a 
production jet engine is explored over a wide range of operating conditions for varying 
rotor imbalance, support stiffness and damping. Interesting dynamical phenomena, 
such as chaos, subharmonic responses, and double-valued responses, are presented and 
discussed. 

KEYWORDS: magnetic bearings, auxiliary bearings, nonlinear rotordynamics. 

NOMENCLATURE 

Cb = auxiliary bearing support damping, ib.s/in. 

F n = normal force, lb 

Fx — external force vector acting on the rotor in X direction 
Fy = external force vector acting on the rotor in Y direction 
I 0 = rotor inertia matrix 

I(b = auxiliary bearing support stiffness, ib/in. 

Kc ~ contact stiffness, lb/in. 

Mb — auxiliary bearing mass, lb.s 2 /in. 

Mk = mass of kth rotor element, lb.s 2 /in. 

N — total number of modes considered 
Nbi = node number at auxiliary bearing #1 
Nb2 = node number at auxiliary bearing #2 
Q,Y = rotor modal coordinate vector in X direction 
Q Y = rotor modal coordinate vector in Y direction 
Rb = radius of auxiliary bearing bore, in. 

Rft = radius of rotor journal, in. 

X* = rotor physical coordinate vector in X direction 
Yft — rotor physical coordinate vector in Y direction 
e = rotor imbalance eccentricity, in. 
g as gravitational acceleration, in./s 2 
t = time, s 

A as deformation at the contact point, in. 

r = ^ r i a ^ 

= rotor free-free modal rotation matrix 
fi = rotor operating speed, rad/s 
^ = rotor free-free modal displacement matrix 
6 = Rb — Rrt) auxiliary bearing clearance, in. 
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f = modal damping coefRcient 

INTRODUCTION 

Active magnetic bearings are one of the most innovative recent developments in the turbomachinery 
field. This technology provides the potential for significant improvements over other types of rotor 
support, including elimination of wear and bearing friction-related energy losses as well as a means 
of actively suppressing rotor vibration. A critical component of any magnetic bearing design is the 
auxiliary bearing, which protects the soft iron core of the magnetic bearing and provides rotor support 
in case of overload or failure of the primary (magnetic) bearing. Magnetic bearing systems appear 
to provide particularly great promise for use in aeronautical applications. In this regard, current 
effort is directed toward developing jet engines with rotors supported by magnetic bearings. For such 
applications, safety is an important concern. Toward this end, it is desirable to design the system to 
operate with auxiliary bearing support for an extended period of time. 

A number of different bearing types have been suggested as auxiliary bearings. These include bushings, 
rolling element bearings, and various types of journal bearings. The most commonly considered are 
rolling element bearings. The major disadvantage associated with using rolling element bearings (or 
bushings) is the requirement of a clearance between the rotor and the inner race of the bearing, without 
which many of the advantages associated with using magnetic bearings would be reduced or eliminated. 
This clearance introduces a nonlinear dynamical feature which may significantly impact the behavior 
of the rotor. 

There are quite a number of studies in the literature concerned with nonlinear rotordynamics. 
Yamamoto (1954) conducted a systematic study of rotor responses involving bearing clearance effects. 
Black (1968) studied the rotor/stator interaction with a clearance. Ehrich (1966, 1988 and 1991), Bently 
(1974), Muszynska (1984) and Childs (1979 and 1982) observed and studied subharmonic responses 
associated with clearance effects. While this work has served to greatly enhance the understanding 
of such systems, more detailed study is needed. Much of this earlier work was conducted from the 
perspective that the clearance is a result of manufacturing error or misfitting and is best eliminated. 
However, in a rotor system fitted with magnetic bearings and auxiliary bearings, the clearance becomes 
a design parameter rather than an irregularity. From tiiis point of view, it is important to develop a 
better understanding of the expected dynamic responses. Such knowledge will provide guidelines for 
the selection of auxiliary bearing parameters. 

There has been relatively little work available in the open literature that is specifically concerned with 
auxiliary bearings. Two papers that are directly related to research on auxiliary bearings are Gelin 
et al., (1990) and Isliii and Kirk (1991). Both of these papers are mainly concerned with transient 
responses. The current work is specifically concerned with developing a better understanding of the 
expected steady-state dynamical behavior for an auxiliary bearing supported rotor. Simulation results 
for a rotor based upon a production jot engine are presented and discussed. 

The model used in the current study has two principal components - the rotor and the auxiliary 
bearings. The rotor is modelled using the free-free bending mode shapes and natural frequencies 
obtained through finite element analysis. The finite element code uses 34 stations and the first four 
modes (two rigid body and two flexible modes) are included in the simulation model. Figure 1 shows 
a schematic diagram of the finite element model. The rotor equations of motion can be expressed in 
terms of modal coordinates as: 


Q A '+2C^nQ.Y + DTQ k +^uQa' 

+ 2fiC^nQy = $ T F.y, (l.a) 

Qy -+2Cw n Qy - flFQ.v +u>;Q,' 

- 2ft<u; B Q.y = $ T F y. 


(1-6) 
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FIGURE 1 DIAGRAM OF THE FINITE FIGURE 2 AUXILIARY BEARING MODEL 

ELEMENT MODEL 


where 


with 


F.v = {FxuF x F Xm }- T , 

Fy = {FViiFyj, — Fy m } _T , 

Q.v = *-‘X R , 

Q y = 

/w ni 0 ... 0 \ 

0 ui n , ... 0 

— 

Vo o ... 

X n ={X M , X m -W -7 \ 

Yr= [Y m ,Y m W T - 

(m = total number of nodes) 


The physical displacements of the rotor at the two auxiliary bearing locations can be obtained using 
the following coordinate transformation: 


N 

XRk = Qxi* 

(k = N BU N B 2) 

y nh = ® ki ® Yi ' 

t=l 


The model for the auxiliary bearings is shown in Figure 2. The governing equations of motion are: 


i^Bk^Bk + Csk^Bk T * A * BkX Bk 

— F nk cos a k + A/flfc g, (2.a) 

MekYek + CsfcVs* + A ' Bk YBk 
~ F nk sina k , 


(2.6) 
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where 


a k = tan 1 


Y Rk - Ysk 

Xnk - Xnk 


(k= N Du N D2 ) 


Fxk = ~ Fnk cos a k + Jl/fc g + Mk e ft 2 cos (Sit), 
Fyk = - F n k a k -f Al k e SI 2 sin (SU). 

The rotor/bcaring interaction is represented with the normal force F nk 

p _ J F c Sk, > 0, 

nk - \ 0, A* < 0, 


where 


Ak - ( XRk - X B k)cosQk + (Vnib - Y B k)sina k - 6 k . 


( 3 ) 


DISCUSSION OF RESULTS 

The model described in the previous section was used to study the steady-state dynamic behavior 
of such a system using direct numerical integration of the equations and a harmonic balance code. 
Studies were performed for varied values of imbalance, support stiffness, and support damping, 
respectively. For the purposes of the current work, the two auxiliary bearings are located at nodes 
3 and 33, respectively. It is assumed that they are identical in terms of stiffness, damping and 
friction characteristics. The nominal system parameters used for the simulation study are Mqi =0.0023, 
Af J 52 :=: 0-0024, K c =2.8550+6, and (=0.03. All the results that are presented correspond to node 3, 
the location of bearing 1. 

The dynamical behavior of a rotor supported by bearings with clearance coupled with the nonsymmetry 
resulting from gravitational effects can be quite complex. The harmonic balance method is first used 
to investigate the synchronous behavior for small imbalance values. (Please note that the complex 
frequency contents associated with medium and large imbalance values makes it a formidable task 
to apply the harmonic balance method for other cases.) Figure 3 is a typical plot of the steady 
state response amplitudes as functions of the rotor speed with and without bearing clearance. The 
nonsymmetric effects resulting from gravity loading is similar to what occurs for a rotor supported by 
nonsymmetric bearing stiffnesses with regard to influence on critical speeds. The first critical speed 
splits into two distinct values. For the X direction, the gravity force tends to keep the rotor in contact 
with the bearing at low operating speed. The apparent stiffness in this direction is approximately the 
same as Kq and the resulting critical speed about- the same as the critical speed for the linear case 
(£=0). For the Y direction, the clearance results in a lower apparent stiffness and a lower critical speed. 
Several higher order pseudo-critical speeds are also created in the operating speed range (about 1500 
rad./scc). It should be noted that the response in the X direction also depart* from the linear case 
at high operating speed. This is because the imbalance force becomes dominant at high rotor speed 
which in turn makes the gravity force less significant and the clearance effect more important. 

Figure 4(a) shows some typical results for varying imbalance. Some of the more dynamically interesting 
results occur for cases of relatively large imbalance. It is observed that imbalance may influence the 
frequency content of the rotor responses quite dramatically at certain operating speeds. There exists as 
many as eight ranges of imbalance values that result in eight different types of rotor responses. In fact, 
subharmonic responses from SI/ 2 through fI/10 are observed. Those subharmonics are not directly 
related to the . y stem’s natural frequencies as were the cases with other researchers' findings (such as 
Ehrich, 1988). Moreover, several types of subharmonic responses may occur for identical parametric 
configurations, but different imbalance values. 

Figures 4(b) and 4(c) show typical results using Kb or Cb as the variable parameter. Clear routes 
to chaos are observed. Kq increases beyond certain value, a period-doubling bifurcation always 
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FIGURE 3 INFLUENCE OF CLEARANCE 
ON CRITICAL SPEEDS 
(e=0.0001 , C B = 150.0, 
K b =0.313xI0 6 ) 



FIGURE 4(a) BIFURCATION DIAGRAM 
WITH IMBALANCE AS THE 
VARIABLE PARAMETER 

(n=1000. <5=0.002, C B = I50, 

K b =0.313x10 6 ) 


takes place. As C b decreases below certain value, the responses always become chaotic. However, the 
bifurcation type is not well defined as for varied Kq- R should be pointed out that even though a lower 
K b may leads to a better system response, it may also fail to protect the magnetic bearings due to the 
fact that it could result in a larger rotor orbit-center offset. In addition, higher bearing damping does 
not necessarily result in synchronous responses. It is observed that, with small imbalance (c < 0.0005), 
the responses are always synchronous for all speeds in the speed range considered (< 1800) if Cq is 
large enough {C b =250.0) and the clearance is small enough (<5 < 0.001). As with the other cases, 
higher bearing damping or lower bearing stiffness or lower imbalance tends to increase the probability 
of synchronous responses. 


.5 


X 




FIGURE 4(b) BIFURCATION DIAGRAM 

WITH BEARING STIFFNESS 

AS THE VARIABLE 

PARAMETER 

(ft=1500, <5=0.002, C B = 150, 

e=0.0008) 


FIGURE 4(c) BIFURCATION DIAGRAM 
WITH BEARING DAMPING 
AS THE VARIABLE 
PARAMETER 
(ft=1500, <5=0.002, e=0.0015, 
K b =0.313x10 6 ) 
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CONCLUSIONS 

As a summary of the results discussed above, the following conclusions can be drawn: 

1. Imbalance may serve to dramatically alter the frequency contents of the rotor responses at certain 
operating speeds. This is particularly evident for cases of large clearance, high bearing stiffness and 
low bearing damping. 

2. For sufficiently high imbalance: 

i. There is a threshold level of damping below which complex dynamical behavior can be expected. 

ii. There is a threshold level of stiffness above which complex dynamical behavior can be expected. 

3. Clear routes to chaos are observed. As bearing stiffness increases beyond a certain value, a period- 
doubling bifurcation takes place which leads to chaos. As bearing damping decreases below a certain 
value, chaos also tends to occur but the bifurcation type is not so clearly defined. 
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ABSTRACT 

Tbis study examines the response of a flexible ro- 
tor supported by load sharing between linear bear- 
ings and an auxiliary clearance bearing. The objec- 
tive of the work is to develop a better understand- 
ing of the dynamical behavior of a magnetic bear- 
ing supported rotor system interacting with auxil- 
iary bearings during a critical operating condition. 
Of particular interest is the effect of coupling be- 
tween the bearing/housing and shaft vibration on 
the rotordynamicaJ responses. A simulation model 
is developed and a number of studies are performed 
for various parametric configurations. An experi- 
mental investigation is also conducted to compare 
and verify the rotordynamic behavior predicted 
by the simulation studies. A strategy for reduc- 
ing sychronous shaft vibration through appropriate 
design of coupled shaft /bearing/ housing vibration 
modes is identified. The results are presented and 
discussed. 

NOMENCLATURE 

C = damping, N-sec/m 

K = stiffness, N /m 

Mb = auxiliary bearing mass, kg. 

Mh = housing mass, kg. 

N — total number of modes considered 

Nb i = node number at leftmost bearing 

Nb 2 = node number at rightmost bearing 

Nbz = node number at auxiliary clearance bearing 

Q s = rotor modal coordinate vector in X direction 

Q y = rotor modal coordinate vector in Y direction 


t = time, s 

X T — rotor physical coordinate vector in X 
direction, m 

Y r — rotor physical coordinate vector in Y 
direction, m 

Xb = auxiliary bearing physical coordinate vector 
in X direction, m 

Yb = auxiliary bearing physical coordinate vector 
in Y direction, m 

Xh = housing physical coordinate vector in X 
direction, m 

Yh. = housing physical coordinate vector in Y 
direction 

A = radial clearance in auxiliary bearing, m 
I a = rotor polar mass inertia matrix, kg-m 2 
'i — rotor free-free modal rotation matrix 
$ = rotor free-free modal displacement matrix 
T = * r I a ¥ 

Q = rotor operating speed, rad/s 
u n = matrix of rotor free-free natural 
frequencies, rad/s 
■tp = imbalance vector 
Subscripts 

bx =■ auxiliary bearing, x-direction 
by = auxiliary bearing, y-direction 
c = contact 

hx = housing, x-direction 
hx = housing, y-direction 
xr = rightmost bearing, x-direction 
xl = leftmost bearing, x-direction 
yr = rightmost bearing, y-direction 
yl = leftmost bearing, y-direction 
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INTRODUCTION 

In recent years, the use of active magnetic bear- 
ings for supporting turbomachinery has been an 
area of interest for both academic researchers and 
turbomachinery users in industry. Magnetic bear- 
ings provide the potential for significant improve- 
ments over other types of rotor support, including 
elimination of wear and bearing friction-related en- 
ergy losses as well as a means of actively suppress- 
ing rotor vibration. However, their use has been 
significantly limited due to a number of technical 
problems. A particular area of concern is the aux- 
iliary bearing, which protects the soft iron core of 
the magnetic bearing and provides rotor support in 
case of overload or failure of the magnetic bearing. 

Typically, the auxiliary bearings have relatively 
small clearances so that magnetic bearing ro- 
tor/stator contact does not occur during bearing 
failure or power loss. Due to these small clear- 
ances, contact between the rotor system and the 
auxiliary bearings can occur during standard oper- 
ation of the magnetic bearings. When this occurs, 
load sharing between the magnetic bearings and the 
auxiliary bearings results, and the rotor system in- 
teracts with its auxiliary bearings. The dynamics 
of such an occurrence must be understood in order 
to properly design the auxiliary bearing system for 
a magnetic bearing supported rotor. 

A number of different bearing types have been 
suggested as auxiliary bearings. These include 
bushings, rolling element bearings, and various 
types of journal bearings. The most commonly 
considered are rolling element bearings. The ma- 
jor disadvantage associated with using rolling ele- 
ment bearings (or bushings) is the requirement of 
a clearance between the rotor and the inner race of 
the bearing, without which many of the advantages 
associated with using magnetic bearings would be 
reduced or eliminated. This clearance introduces 
a nonlinear dynamical feature which may signifi- 
cantly impact the behavior of the rotor. 

There are quite a number of studies in the lit- 
erature concerned with nonlinear rotordynamics. 
Ehrich (1965 and 1967) studied the rotor/stator 
interaction with a clearance and found zones of 
bistable synchronous behavior. Black (1968) also 
found these bistable interaction zones. Ehrich ex- 
tended this research by predicting subharmonic 
behavior (1966 and 1988), superharmonic behav- 
ior (1992), and chaotic responses (1991). Bently 
(1974), Muszynska (1984), and Childs (1979 and 
1982) also observed and studied subharmonic re- 
sponses due to clearance effects. 

There is relatively little work available in the 
open literature that is specifically concerned with 
the dynamics between the rotor system and the 
auxiliary bearings. The research that has been 


performed to date is primarily concerned with the 
transient behavior of rotor drop on the auxiliary 
bearings due to power failure or inoperative mag- 
netic bearings. Two papers with this focus are 
Gelin et al., (1990) and Ishii and Kirk (1991). 

The present work is concerned with developing 
an understanding of the dynamic behavior of a 
rotor system supported by load sharing between 
the linear bearings and an auxiliary bearing with 
clearance. Of particular interest is the influence of 
coupled shaft /bearing/housing vibration modes on 
the rotordynamical behavior of such a system. 

SIMULATION MODEL 

Figure 1 shows the simulation model used for 
this investigation. It consists of a flexible rotor 
supported at both ends by magnetic bearings. A 
rigid disk with adjustable imbalance is placed at the 
midpoint of the bearing span. As a simplification, 
the magnetic bearings are modeled as spring and 
dashpot systems and interaction with a single 
auxiliary bearing is considered. The auxiliary 
bearing was modeled as an antifriction bearing with 
a clearance and a mass attached to a housing. The 
housing in turn has a mass, stiffness, and damping. 
Figure 2 shows the auxiliary bearing model used in 
the simulation. 


Rigid Disk 



Figure 1 Schematic Diagram of Simulation Model 

The rotor is modelled using the free-free bend- 
ing mode shapes and natural frequencies obtained 
through finite element analysis. The finite element 
code uses 19 stations and the first four modes (two 
rigid body and two flexible modes) are included in 
the simulation model. 

Using these simplifications, the equations of mo- 
tion for the system can be written as follows: 

Qx+flFQy + + ^^'(F i'lin + F c,nl) 


2 



= * r F x , imiM (l.«) 

Qy — firQ r + U)^Qy 4" ^ (-fj y,lin 4" Fy,nl} 



— 4? F y 

(1.6) 

M b X b 4 

-^(X b - Xh.) 4- K bx (X b - X h ) 



= F r,ni* 

(l.c) 

M b Y b 4 

~(Y b - Yfc) + K by (Y b - Y fc ) 



= Fy jn i , 

(l.d) 


M h X h +(^ + ~~)Xh 4- (iTox + K hx )X h 

= C&x-Xfc + KbxX 6, (l* e ) 

M fc y fc 4-(^ 4- ^)Y> + (/V6v + 

+ (i./) 


where 

F x ,«„ = * Nil (K xl Q x 4- ^Qr) 

+ * W - w (lT, r Q, + ^Q«) 

Fy T jj n = ^jVii(-^ylQy H q” Q y) 

+ $Nt 3 {KyrQy + -^pQy) 

= i>n 2 cosnt, 

Fy f ijnb = sinCtt^ 

F^ = ^-A)( $Q y^ ) 

F y , n ( = - A)( #Q y~ — ) 

6 = ^/(*Qx - Xb) 2 + (*Qy - n) 2 , 
<f> = 1 i/ 6 > A 
0 otherwise 

Q* = * -1 X r , 

Qy = ^ _1 Y r , 

with 

X r = {X rl ,X r2 ,...,X rm } T , 

Y r = {y rl ,y r 2 ,...,y m } r . 

(m = total number of nodes) 


The physical displacements of the rotor at the 
auxiliary bearing or magnetic bearing locations 
can be obtained using the following coordinate 
transformation: 

jv 

Xrk = ^ ^ ^ ki Q xi » 

*“ l (k = N bu N b2 ,N b 3 ) 

y-fc = ^ ^ fci Q yt > 

t=l 


Rotor 



Figure 2 Auxiliary Bearing Model 


EXPERIMENTAL MODEL 

Experimental tests were performed in order to 
validate the behavior predicted by the simulation 
model and to gain some insight into the dynamical 
responses that are to be expected. A schematic 
diagram of the rotor test rig that was used in the 
experimental work is shown in Figure 3. 

The rotor used in this study has two basic 
components: a flexible shaft and an auxiliary 

clearance bearing. The shaft is made of steel 
and is 0.374 inches in diameter and 18.0 inches 
in length. It is supported at 1.0 inch from the 
right end by ball bearings suspended in a frame 
by four springs and at 1.0 inch from the left end by 
a bushing with a tight clearance. These supports 
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Figure 3 Experimental Model 


represent the magnetic bearings. The stiffness of 
the left support is 17511 N/m (101.5 lb/in) for 
both the horizontal and vertical directions. It 
is used to somewhat isolate the rotor from the 
effects of the flexible coupling which attaches the 
rotor to the motor and to enforce low amplitude 
vibration at this location to protect the motor. The 
stiffness of the right support, for both horizontal 
and vertical directions, is 2539 N/m (14.5 lb/in). 
This lower stiffness allows for significant vibration 
of the rotor in the speed range of the motor. A 
rigid disk with an adjustable imbalance is placed 
at the midpoint of the bearing span. The auxiliary 
clearance bearing/housing consists of a bushing 
suspended in a frame by four springs (Figure 4). 
It is situated at the right end of the rotor. The 
clearance is adjustable by changing the bushing. 
The auxiliary bearing/housing stiffnesses are varied 
by interchanging the springs. 

The rotor is driven by an adjustable speed mo- 
tor with feedback speed controller. Shaft vibration 
is measured using eddy current proximity displace- 
ment sensors fixed so as to measure displacement 
in the vertical and horizontal directions. The dis- 
placement signals were recorded and analyzed with 
a signal analyzer. 
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DISCUSSION 

A series of parallel studies were conducted using 
the simulation and experimental models discussed 
in the preceding sections. The following discussion 
is aimed at comparing the predictions of the simula- 
tion model with experimentally observed responses, 
with the objective of obtaining insight into the be- 
havior of flexible rotor system dynamics due to load 
sharing between the magnetic bearings and an aux- 
iliary bearing. 

Using the data from measurements and calcula- 
tions, the stiffness, damping, and mass character- 
istics of the experimental rotor system were identi- 
fied. The imbalance configuration consists of 3.0 x 
10~ 6 kg-m located at the left end of the rotor shaft 
(from the flexible coupling to the electric motor) 
and a variable imbalance located on the rigid disk. 
Linear, hysteretic, and coulomb friction damping 
models were considered. The support damping ap- 
pears to result principally from hysteresis. If syn- 
chronous vibration is assumed, this results in a 
damping coefficient scaled by the rotor speed. The 
numerical values are given in Table 1. 

The governing parameters for an auxiliary bear- 
ing are stiffness, mass, damping, clearance, and ax- 
ial location. The axial location generally is fixed 
by other considerations, such as space requirements 
and the need to be close to the magnetic bearing 
to better protect it. This study focusses on the 
remaining design parameters. For the system stud- 
ied, significant nonsynchronous vibration occured 
only for the 2u> component and for a limited ro- 
tor speed range. Figure 5 illustrates the observed 
behavior for various imbalance values. The ampli- 
tudes are quite small and, interestingly, virtually 
independent of imbalance. The remainder of this 
discussion will focus on the synchronous response 
amplitudes. 
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Parameter | Value | 

Units 

K b 

471 

N/m 

K h 

1.8x10' 

N/m 

K c 

87.557 

N/m 

K !x , K)y 

17,510 

N/m 

K«, K ry 

2,539 

N/m 

Cjx? Cjy 

1.000 

N/m 

CfX? C r y 

500 

N/m 

Cly> C|y 

1,000 

N/m 

Cry ? Cry 

500 

N/m 

Cbx, Cby 

2,000 

N/m 

Chx» Chy 

0 

N/m 

Mb 

0.033 

kg 

M h 

0.296 

kg 

i>\ 

j 3.0 x 10~ 6 

kg-m 

i>\o 

2.5 x 10-* 

kg-m 

A 

5.0xl0 -i 

m 


Table 1 Simulation Model Parameters 



Figure 5 Experimental Results for 
2u Component 

The simulation responses discussed below were 
determined from a single-term harmonic balance 
analysis that was numerically implemented and ver- 
ified at selected points using direct numerical inte- 
gration of the governing equations of motion. The 
harmonic balance procedure is described by Nayfeh 
and Mook (1979). Figures 6. a - 6.e show the re- 
sponse amplitudes predicted from the simulation 
study for the parameters of Table 1. Figures 7.a - 
7.c show some corresponding experimental results. 
The rotor speed range extends to above the second 
critical speed, with the first critical speed primar- 


ily a rotational mode and the second critical speed 
primarily a translational mode. 

Figures 6. a presents simulation results for vary- 
ing auxiliary bearing stiffness and Figure 7.a shows 
the corresponding experimental values. There is 
relatively good agreement between the actual and 
predicted responses except for the peak values 
which are very sensitive to damping. At lower ro- 
tor speeds, the stiffness effect from the auxiliary 
bearing tends to dominate and serves to increase 
the effective lower critical speed. At higher ro- 
tor speeds, the inertial effects of the bearing tend 
to dominate and the added mass serves to lower 
the effective second critical speed. This effect is 
balanced somewhat for higher stiffnesses (as is ex- 
pected), as can be seen for the case with Kb=2,625 
X /m. The damping added to the system from the 
auxiliary bearing serves to attenuate the peak re- 
sponse amplitudes for both critical speeds. The re- 
sponse amplitudes between the two response peaks 
are affected little by the auxiliary bearing, with the 
rotor and bearing/' housing vibrations tending to 
decouple for response amplitudes below the clear- 
ance. Above the second critical speed, the pres- 
ence of the auxiliary bearing serves to reduce the 
response amplitudes for the remainder of the plot. 
For this region, there is significant coupling be- 
tween the rotor/bearing/housing vibration and the 
response amplitudes are below the clearance for 
speeds greater than about 550 rad/sec. The ro- 
tor speed ranges for which coupling of the ro- 
tor/bearing/housing vibration will occur depend 
upon the phase of the combined system response, 
as described by Black (1968). For the purposes of 
this study, Kb=471 N/m is chosen as a reasonable 
value and variations of other parameters are per- 
formed with this auxiliary bearing stiffness. 



Figure 6. a Simulation Results for Varying 
Auxiliary Bearing Stiffness 
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Figure 6.b Simulation Results for Varying 
Auxiliary Bearing Mass 



w (rad/s) 

Figure 6.c Simulation Results for Varying 
Auxiliary Bearing Damping 

Figures 6.b shows the predicted rotor response 
for the parameters of Table 1 with the bearing mass 
varied. These results show that the more massive 
bearing/housing configurations can result in much 
lower amplitude rotor responses. In fact, for suffi- 
ciently high mass values, the response amplitude is 
dramatically reduced for the rotor speed range from 
about 600 rad/sec to about 780 rad/sec where it 
merges with other amplitude response curves. Fig- 
ure 6.c shows the effect of auxiliary bearing damp- 
ing on the rotor response. While the response am- 
plitudes remain relatively unchanged for the ma- 
jority of rotor speeds, higher damping does serve 
to significantly reduce the rotor response ampli- 
tudes through the critical speeds. Damping also 
effects the phase characteristics of the combined 
rotor /bearing/housing system and changes the ro- 



w (rad/s) 

Figure 6.d Simulation Results for Varying 
Imbalance 



Figure 6.e Simulation Results for Varying 
Auxiliary Bearing Clearance 

tor speed ranges for which coupled responses occur. 
Figures 6.d and 7.b show the rotor response for the 
parameters of Table 1 with the imbalance varied. 
As expected, the response amplitudes increase with 
increasing imbalance. Figures 6.e and 7.c show the 
rotor response amplitudes with the auxiliary bear- 
ing clearance varied. As the clearance increases, the 
dynamics of the combined rotor /bearing/housing 
vibration tend to decouple over a wider speed range 
between the two critical speeds. Again, such be- 
havior is expected from the analysis work for a ro- 
tor operating in an annular clearance performed by 
Black (1968). 

Close examination of the above results reveals a 
very interesting trend. Turbomachine are designed 
to operate at rotor speeds between the various sys- 
tem critical speeds so as to minimize synchronous 
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Figure 7.a Experimental Results for Varying 
Auxiliary Bearing Stiffness 



Figure 7.b Experimental Results for Varying 
Imbalance 

response amplitudes. For the present system, the 
lowest amplitude responses (for operation above 
the second critical speed) occur for the auxiliary 
bearing configurations where there is the largest de- 
crease of the second effective critical speed, which 
is due to a coupling of rotor/stator/housing vibra- 
tion. While not really suprising, this result indi- 
cates an interesting strategy for the design of aux- 
iliary support using clearance bearings. Selection 
of combined shaft/auxiliary bearing/housing vibra- 
tion characteristics such that the effective critical 
speeds nearest the operating speed (directly above 
and below) shift away from it should result in lower 
amplitude rotor vibration and provide the best op- 
portunity to protect the magnetic bearing. How- 
ever, the phase characteristics of the combined sys- 
tem must be such that vibration coupling occurs 


Figure 7.c Experimental Results for Varying 
Auxiliary Bearing Clearance' 

and care must be taken to ensure that this is the 
case if the maximum benefit is to be derived from 
the auxiliary bearing. 

CONCLUSIONS 

A study of the dynamical behavior of a flexible 
rotor supported by linear bearings (representing a 
set of magnetic bearings) and an auxiliary bearing 
with clearance has been presented. Parallel sim- 
ulation and experimental studies have been per- 
formed for a variety of parametric configurations. 
The influence of bearing mass, stiffness, damping, 
clearance, and imbalance were examined. It was 
observed that the dynamic behavior of a rotor in- 
teracting with an auxiliary bearing depends very 
strongly on the structural parameters of the auxil- 
iary bearing and associated housing. Appropriate 
selection of these parameters is critical if acceptable 
vibration characteristics are to be obtained for such 
systems. In general, one must carefully consider the 
influence of housing/bearing dynamics and how to 
best take advantage of favorable coupled modes of 
vibration in the development of auxiliary bearing 
designs. 
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ABSTRACT 

This paper investigates the steady-state responses of a 
rotor system supported by auxiliary bearings in which there 
is a clearance between the rotor and the inner race of 
the bearing. A simulation model based upon the rotor 
of a production jet engine is developed and its steady- 
state behavior is explored over a wide range of operating 
conditions for various parametric configurations. Specifically, 
the influence of rotor imbalance, clearance, support stiffness 
and damping is studied. Bifurcation diagrams are used as 
a tool to examine the dynamic behavior of this system as a 
function of the afore mentioned parameters. The harmonic 
balance method is also employed for synchronous response 
cases. The observed dynamical responses is discussed and 
some insights into the behavior of such systems are presented. 

NOMENCLATURE 

Cb = auxiliary bearing support damping, lb.s 2 /in. 

F n = normal force, lb 

Fjf = external force vector acting on the rotor in X direction 
Fy = external force vector acting on the rotor in Y direction 
I a = rotor inertia matrix 

Kb — auxiliary bearing support stiffness, lb/in. 

Kc = contact stiffness, lb/in. 


Mb — auxiliary bearing mass, Ib.s 2 /in. 

Mk — mass of kth rotor element, lb.s 2 /in. 

N = total number of modes considered 
NB 1 = node number at auxiliary bearing #1 
NB2 = node number at auxiliary bearing #2 
Qx — rotor modal coordinate vector in X direction 
Qy = rotor modal coordinate vector in Y direction 
Rb = radius of auxiliary bearing bore, in. 

R m = radius of auxiliary bearing pitch, in. 

Rr ss radius of rotor journal, in. 

X.r = rotor physical coordinate vector in X direction 
Yr = rotor physical coordinate vector in Y direction 
e = rotor imbalance eccentricity, in. 
g = gravitational acceleration, in./s 2 
t = time, s 

A = deformation at the contact point, in. 

r = * T i a * 

* = rotor free-free modal rotation matrix 
fl = rotor operating speed, rad/s 
u> n = rotor natural frequency, rad/s 

= rotor free-free modal displacement matrix 
S = Rr — Rr, auxiliary bearing clearance, in. 
fi = dynamic friction coefficient 
= rolling friction coefficient 

ipB = angular displacement of auxiliary bearing inner-race 
£ = modal damping coefficient 
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INTRODUCTION 

One of the most innovative developments in the turbo- 
machinery field involves the use of active magnetic bearings 
(AMB) for rotor support. This technology provides the po- 
tential for significant improvements in the dynamic behavior 
of rotor systems, allowing for loading, eccentricity, shaft po- 
sition and vibration to be continuously monitored and con- 
trolled. In order to protect the soft iron cores of the magnetic 
bearings and to provide rotor support in the event of failure 
of the bearing or during an overload situation, auxiliary (or 
backup) bearings are usually included in the rotor design. 
There are a number of auxiliary bearing configurations that 
have been considered. Typically a clearance bearing (bushing 
or rolling element bearing), in which there is a gap between 
the rotor and the inner race of the bearing, is used. The clear- 
ance introduces a nonlinear dynamical feature which may 
significantly impact the behavior of the rotor. 

There are a number of studies in the literature concerned 
with the dynamics of rotors with clearance effects. Ya- 
mamoto (1954) conducted a systematic study of rotor re- 
sponses involving bearing clearance effects. Black (1968) 
studied the rotor/stator interaction with a clearance. He 
concluded that rotor /stator interactions may occur in a vari- 
ety of forms and circumstances, including jump phenomena. 
Ehrich (1966) reported the first identification of a second 
order subharmonic vibration phenomenon in a rotor system 
associated with bearing clearance (1966). Bently (1974) pub- 
lished experimental observations of second and third order 
subharmonic vibration in a rotor system. Later, Muszynska 
(1984) cited the occurrence of second, third, and fourth order 
subharmonic responses in a rotor rubbing case and Ehrich 
(1988 and 1991) observed eighth and ninth order subhar- 
monic vibration as well as chaotic vibration in a high speed 
turbomachine. Childs (1979 and 1982) published two papers 
to explain the mechanism for the second and third order 
subharmonic responses noted above. Choi and Noah (1994) 
investigated the influence of bearing clearances on rotordy- 
namic responses, including bifurcation behavior and chaos. 
This paper also contains an excellent review and list of ad- 
ditional references on the topic. 

There has been quite a bit of recent work that is specifi- 
cally concerned with auxiliary bearings in magnetic bearing 
supported rotor systems. Gelin et al. (1990) studied the 
transient dynamic behavior of rotors on auxiliary bearings 
during the coast down. Ishii and Kirk (1991) and Kirk and 
Ishii (1993) investigated the transient responses of a flexi- 
ble rotor during the rotor drop after the magnetic bearings 
become inactive. Schweitzer, et al. (1994) presents a good 
discussion of issues related to the touch-down dynamics of 
rotors on auxiliary bearings. While these studies have greatly 
enhanced the understanding of the dynamics of rotors sup- 
ported by auxiliary clearance bearings, most of them have 


been performed from the perspective that the rotor will be 
shut down if failure of one or more of the magnetic bearings 
occurs. As a result, work in this area has concentrated on 
transient dynamical behavior for a rotor with a decreasing 
spin speed. From a practical perspective, such an assump- 
tion is appropriate for noncritical applications, such as power 
generators and compressors. However, magnetic bearing sys- 
tems appear to provide particularly great promise for use in 
aeronautical and aerospace applications. There are active 
programs at many of the major jet engine manufacturers to 
develop engines supported by magnetic bearings. Safety is a 
major concerns in any aeronautical design. In particular, it 
is desirable for an aircraft to be able to return to its base and 
land if the magnetic bearing(s) fail. In addition, flywheel en- 
ergy storage systems which require extended operation of the 
rotor on auxiliary bearings are being considered for space ap- 
plications. Toward this end, it is desirable to design the rotor 
system to take maximum advantage of the backup bearings 
and use them as true auxiliary bearings to provide support 
during critical situations in a safe and consistent manner. If 
clearance bearings are to be used for such applications, an 
important concern in this regard is the dynamic behavior of 
the rotor when it is in contact with the auxiliary bearing. 
If safe and effective operation is to be ensured during these 
periods, it is essential that designers have a very good under- 
standing of the steady-state dynamics of rotor systems with 
clearance effects. 

In the current paper, simulation studies are performed for 
a complex rotor system supported by auxiliary bearings with 
clearance at each end of the rotor. This work is specifically 
concerned with systems in which the clearances are quite 
small (on the order of a few mils), which is appropriate 
for applications in which the backup bearing is acting to 
provide rotor support on a consistent basis. The influence 
of rotor imbalance, clearance, support stiffness and support 
damping are investigated using direct numerical integration 
of the governing equations of motion and the harmonic 
balance method. Some insights are obtained with regard 
to the frequency and amplitude behavior of the steady-state 
vibration of such a system. Based upon these results, some 
design guidelines and suggestions are presented. 

SIMULATION MODEL 

The rotor is modelled using free-free normal mode shapes 
and natural frequencies obtained through finite element 
analysis. The model data is representative for the rotor 
of a jet engine. Fig. 1 shows a schematic diagram of the 
FEM rotor model. The torsional motion of the shaft is not 
considered in this paper. Using state space representation 
and modal coordinates, the equations of motion for the rotor 
are expressed as 
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The rotor/bearing interaction is represented with the normal 
force Fnk 


Qjr+2Cw„Qjr + firQy 

+ 2QC w nQy = 3> T F.y, (l.o) 

Qy+2Cw ra Qy — ^I?Qx +w^Qy 

-2nCw„Qx = $ T Fy, (1.6) 

where 

Fx = {Fxi,Fx2, • • • , Fy m } -1 , 

Fy = {Fyi, Fyi, . . . , Fym } -1 , 

Qx = Q-'Xr, 

Qy = $~ 1 Yr, 

with 

(m = total number of nodes) 


The physical displacements of the rotor at the two aux- 
iliary bearing locations can be obtained using the following 
coordinate transformation: 


N 

XrIc = $ik Qxi i 

i=l 

N 

Yric = y^tfc Qn, 

»= i 


(k — NB1,NB2) 


The equations of motion for the auxiliary bearings are 
derived using the model shown in Fig. 2 


MBkXsk + CRkXsk + KRkXsk 

= Fnkcosatk — F t k sinak + MBkg, (2.a) 
MBkYBk + CBkYBk + KskYsk 

= F nk sin a k + F ik cos a k , (2.6) 

JBki>Bk + C Brjik^Bk 

= FtkRsk — ftipFnkRmk, (2.c) 


where 


ak = tan 1 


Ym-Ysk 
XRk — Xsk 


(k = NBl,NB2) 


The forces resulting from gravity, imbalance, and ro- 
tor/bearing contact are expressed as follows. 


Fxk = - F n k cos ak + Ftk sin a h 
+ Mk g + Mk e fi 2 cos (fit), 
Fvh = - Fnk sin a k - F tk cos a k 
+ M k e£2 2 sin 


where 


F n k = 


f K c S k , 

\0, 


Afc < 0, 
Afc >0, 


(3-a) 


A * = (X R k — Xsk ) cos + (Y Rk — Ysk) sin at — 6* 

and the Coulomb friction force Ftk • As long as there exists 
slip at the contact point, the friction force obeys 

Ftk = p F n k- (3.6) 

However, when there is no slip at the contact point, the 
friction forces are solved from equations (1) and (2) using 
the kinematic constraint that the circumferential velocities 
of the rotor and the inner-race of the back-up bearing at 
the contact point equal to each other. At the same time, if 
this solved friction force exceeds the maximum static friction 
force (= F n k) j equation (3.b) applys again. 


DISCUSSION OF RESULTS 

The rotor is modeled with 34 stations (as shown in 
Figure 1) and the first four modes (two rigid body and 
two flexible modes) are included in the simulation model. 
The two auxiliary bearings are located at nodes 3 and 33, 
respectively. This arrangement is taken to represent one 
of the most technically feasible configurations in that it 
greatly simplifies bearing maintenance. It is assumed that 
the two auxiliary bearings are identical in terms of stiffness, 
damping and friction characteristics. Some nominal system 
parameters used for the simulation study are Mr =0.0023, 
i^(7=2.855e-b6, and £=0.03. To avoid excessive cluttering 
of the plots, all the results that are presented in this paper 
correspond to the location of the first bearing, which is node 
number 3. 

The primary design parameters for an auxiliary clearance 
bearing are stiffness, damping, clearance, and load level (im- 
balance). It is fairly obvious that rotor responses involv- 
ing nonsymmetric bearing clearance effects are very complex 
problems and it is quit difficult to obtain a global picture of 
the system responses with numerical integration alone. In the 
first phase of the present work, the harmonic balance method 
is used to investigate the global system behavior. Please 
note that the complex frequency contents associated with 
medium and large imbalance values makes it a formidable 
task to apply the harmonic balance method for configura- 
tions other than for purely synchronous behavior. However, 
an adequately balanced rotor system should have very small 
imbalance under normal conditions and, thus, should exhibit 
primarily synchronous responses. 
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Variation of rotor response amplitudes with rotor speed 
for the system with and without bearing clearance is shown 
in Figure 3(a). Examination of the response curves indicates 
that the effect of nonsymmetries resulting from gravitational 
loading, coupled with bearing clearance, is similar to that 
of asymmetric support stiffnesses with regards to critical 
speeds. The first critical speed actually splits into two 
pseudo-critical speeds. In the X direction, the gravity force 
tends to keep the rotor in contact with the bearing at low 
operating speed. Thus, the apparent stiffness is almost the 
same as Kb and the pseudo-critical speed is nearly the same 
as the critical speed for the linear case (8 = 0). In the Y 
direction, the clearance results in a lower apparent stiffness 
and, consequently, an additional lower-value pseudo-critical 
speed. It is seen that several higher order additional pseudo- 
critical speeds are created in the operating speed range in 
addition to the 1st additional pseudo-critical speed. The 
response in the X direction also departs from the linear case 
at high operating speed. This is because the imbalance force 
becomes dominant at high rotor speed which in turn makes 
the gravity force less significant and the clearance effect more 
important. An interesting phenomena that can result as 
the clearance is increased is that of double-valued responses. 
Figure 3(b) shows double-valued responses in the Y direction 
for four different values of clearance. It is seen that a larger 
clearance results in wider rotor speed range of double- valued 
responses. It is also observed that as clearance increases, 
the apparent stiffness decreases and the first pseudo-critical 
speed shifts to a lower value. Figure 3(c) shows the double- 
valued responses in the X direction. Even though the jumps 
themselves are smaller in magnitude, they are more obvious 
in trend. Notice how little the change is for the first pseudo- 
critical speed in the X direction. The system behaviors for 
higher operating speed range are not shown in Figures 3.b 
and 3.c so that the jump phenomena can be more clearly 
illustrated. The system responses in the high operating speed 
range with the same parameters exhibit primarily amplitude 
changes. 

In the second phase of this study, bifurcation diagrams 
are constructed using the steady-state system responses ob- 
tained through direct numerical integration of the govern- 
ing equations of motion. Such an approach allows for the 
study of medium and high imbalance cases for which non- 
synchronous responses may occur. Near-zero initial condi- 
tions were used, simulating situations where the AMBs are 
functioning properly prior to a system failure. Bifurcation 
diagrams for varying bearing stiffness, bearing damping, and 
clearance are shown. 

Bifurcation diagrams with bearing stiffness as the variable 
parameter are shown in Figure 4 and bifurcation diagrams 
with bearing damping as the variable parameter are shown 
in Figure 5. These response characteristics are typical for 


a quite wide range of configurations. There is a threshold 
stiffness value (about 2.5 x 10 5 for the particular case 
shown) above which the behavior is always multi-frequency. 
Conversely, there is a threshold damping value (about 150 for 
the particular case shown) below which the behavior tends 
to always be multi-frequency. Figure 6 shows bifurcation 
diagrams with clearance as the variable parameter. Two 
values of imbalance are shown. For the lower value (Figure 
6. a and Figure 6.b), the responses are synchronous over a 
wide range of clearance. However, there is a region (0.002< 8 
0.0032) for which the behavior is chaotic. However, for the 
larger value of imbalance (Figure 6.c and Figure 6.d) it is 
observed that the there is a threshold value of imbalance (6 > 
0.0024) for which the behavior is always chaotic. Studies 
were also conducted using rotor imbalance as the variable 
parameter. The variation in response characteristics for this 
case is quite similar to that for the case of varied clearance 
and are not shown for the purpose of brevity. 

A common feature among all of the observed responses is 
that for very small imbalance, the responses are always syn- 
chronous. The imbalance range that result in synchronous 
dominated responses depends on several system parameters. 
For low back-up bearing stiffness (such as Kb = 0.213e+6) 
and moderate damping ( Cb = 157 . 0 ), the responses are 
almost always synchronous. Only ft/2 subharmonic are ob- 
served at a few operating speeds with a very narrow range of 
imbalance. It should be noted that even though a lower Kb 
may leads to a better system response, it may also fail to pro- 
tect the magnetic bearings due to the fact that it could result 
in a larger rotor orbit-center offset. However, the dramatic 
response changes discussed above may occur again if the 
damping becomes small (such as Cb = 50 . 0 ) even though 
the stiffness still remains small. On the other hand, increas- 
ing the damping Cb alone may not be able to eliminate those 
dramatic changes. It is observed that those changes can still 
occur for Cb being as large as 700.0. Reducing the size of 
clearance 8 may not eliminate the multi-frequency behavior 
at certain speeds, but it does tend to narrow the operat- 
ing speed range where those changes occur. For example, 
multi-frequency behavior is eliminated for ft > 1500 when 
8 is reduced from 0.002 to 0.001 with all other parameters 
remaining the same, but such behavior may still occur for 
ft < 1400 . 

CONCLUSIONS 

The harmonic balance method and direct numerical inte- 
gration has been used to study the response characteristics 
of a jet engine rotor supported by bearings with a clearance. 
Variations of bearing stiffness, bearing damping , clearance, 
and rotor imbalance were considered. The dynamic reponses 
are described in some detail and some suggestions are pro- 
vided with regard to the design of an appropriate auxiliary 
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bearing configuration for this rotor system. The basic ob- 
servations are generally applicable to any rotor system. In 
summary, the following observations have been made. 

1. A close examination of the rotor behavior reveals three 
basic response modes. The first is primarily synchronous 
oscillation at the bottom of the bearing dominated by 
gravitional loading holding the rotor against the bearing. 
As vibration amplitudes increase, the motion transitions to 
a a bouncing mode with multiple frequency components. 
Finally as the vibration amplitudes become sufficiently 
large, the motion transitions back to purely synchronous 
whirling. 

2. Steady-state synchronous, multi-synchronous, sub- syn- 
chronous, and chaotic responses may occur and are ob- 
served. Gravitational effects produce a side force on the 
rotor (a nonsymmetric effect) that may be beneficial (by 
keeping the rotor in contact with the bearing inner race) or 
may be harmful (by encouraging the development of multi- 
frequency vibrations /chaos) depending upon the particu- 
lar system configuration. Such effects are certainly present 
in real rotors and need careful consideration in design anal- 
yses for such systems. 

3. Imbalance may serve to dramatically alter the frequency 
contents of the rotor responses at certain operating speeds. 
This is particularly evident for cases of large clearance, 
high bearing stiffness and low bearing damping. Generally, 
auxiliary bearing designs with low clearance, low support 
stiffnesses, and high support damping tend to discourage 
the development of multi-frequency /chaotic behavior and 
provide the most favorable rotordynamic behavior. How- 
ever other requirements, most particularly the requirement 
that no contact occur between the soft iron components 
of the AMB, may impose additional constraints that must 
be considered. 
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Figure 1: Diagram of Finite Element Rotor Model 
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Figure 2: Auxiliary Bearing Model 
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Figure 3.a: Synchronous Response Amplitudes with 
Auxiliary Bearing Clearance as the 
the Variable Parameter 
C B = 150, e = 0.0001, K b = 0.313e6 
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Figure 3.b: Synchronous Response Amplitudes with 
Auxiliary Bearing Clearance as the 
the Variable Parameter 
C B = 300, e = 0.0001, K B = 0.313e6 
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Figure 3.c: Synchronous Response Amplitudes with 
Auxiliary Bearing Clearance as the 
the Variable Parameter 
C B = 300, e = 0.0001, K B = 0.313e6 
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Figure 4.a: Bifurcation Diagram with Auxiliary Bearing 
Stiffness as the Variable Parameter 
G b = 250, cj = 1500, e = 0.0008, d = 0.002 
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Figure 4*b: Bifurcation Diagram with Auxiliary Bearing 
Stiffness as the Variable Parameter 
C B = 250, u = 1500, e = 0.0008 d = 0.002 
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Figure 4.c: Bifurcation Diagram with Auxiliary Bearing 
Stiffness as the Variable Parameter 
C B = 250, u = 1500, e = 0.0004 d = 0.002 
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Figure 4.d: Bifurcation Diagram with Auxiliary Bearing 
Stiffness as the Variable Parameter 
C B = 250, lj = 1500, e = 0.0004 d = 0.002 
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Figure 5. a: Bifurcation Diagram with Auxiliary Bearing 
Damping as the Variable Parameter 
K b = 0.313e6, u = 1500, e = 0.0015 d = 0.002 
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Figure 5.b: Bifurcation Diagram with Auxiliary Bearing 
Damping as the Variable Parameter 
K b = 0.313e6, u = 1500, e = 0.0015 d = 0.002 
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Figure 5.c: Bifurcation Diagram with Auxiliary Bearing 
Damping as the Variable Parameter 
K b = 0.313e6, u = 1000, e = 0.0015 d = 0.001 
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Figure 5.d: Bifurcation Diagram with Auxiliary Bearing 
Damping as the Variable Parameter 
K b = 0.313e6, cj = 1000, e = 0.0015 d = 0.001 
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Figure 6.b: Bifurcation Diagram with Auxiliary Bearing 
Clearance as the Variable Parameter 
C B = 250, K b = 0.313e6, w = 1500, e = 0.0005 
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Figure 6.c: Bifurcation Diagram with Auxiliary Bearing 
Clearance as the Variable Parameter 
C B = 250, K b = 0.313e6, w = 1500, e = 0.0010 
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Figure 6.d: Bifurcation Diagram with Auxiliary Bearing 
Clearance as the Variable Parameter 
C B = 250, K b = 0.313e6, t o = 1500, e = 0.0010 
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- ABSTRACT 

Auxiliary bearings are a critical feature of any magnetic bearing system. They protect 
the soft iron core of the magnetic bearing during an overload or failure. An auxiliary 
bearing typically consists of a rolling element bearing or bushing with a clearance gap 
between the rotor and the inner race of the support. The dynamics of such systems can be 
quite complex. It is desired to develop a rotordynamic model which describes the dynamic 
behavior of a flexible rotor system with magnetic bearings including auxiliary bearings. The 
model is based upon an experimental test facility. Some simulation studies are presented 
to illustrate the behavior of the model. In particular, the effects of introducing sideloading 
from the magnetic bearing when one coil fails is studied. These results are presented and 
discussed. 

NOMENCLATURE 
C = damping, N-s/m 

Cl = clearance of auxiliary bearing, m 
D = nominal gap thickness, m 
Fx — external force vector acting on the rotor in 
X direction, N 

F y = external force vector acting on the rotor in 
Y direction, N 

G = gravitational acceleration, m/s 2 

I = rotor inertia matrix 

i = current, amp 

K = stiffness, N/m 

k = gain value 

L = equivalent circuit length 

M = mass, kg 

m = total number of nodes 

N = total number of modes considered 

Nbi = node number at rotor left end 

N B 2 — node number at AMB 

Nb 3 = node number at auxiliary bearing 
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Ni m b = node number at imbalance location 
Nt = number of turns of wire per coil 
Qx — rotor modal coordinate vector in X direction 
Qv = rotor modal coordinate vector in Y direction 
Rr = radius of rotor journal, m. 

Rs = radius of auxiliary bearing bore, m. 
t = time, s 
tf = fail time, s 
v = voltage, volt 

v r = relative velocity at auxiliary bearing 
rotor /stator contact point 
X# = rotor physical coordinate vector in 
X direction 

Xs = stator physical coordinate vector in 
X direction 

= rotor physical coordinate vector in 

Y direction 

Y 5 = stator physical coordinate vector in 

Y direction 

a = acceleration, rad/sec 2 
6 = dynamic clearance 

r = 

— rotor free- free modal displacement matrix 
\l> = rotor free-free modal rotation matrix 
ij> = imbalance vector 

= permeability of free space 
ji r = permeability of silicon steel 
= rotor operating speed, rad/s 
r = time constant 
9 — angular position of the shaft 
£ = modal damping coefficient 

Subscripts 
a = first 

5 = auxiliary bearing 
b = second 
C = contact 
= gain constant 
/i = horizontal 
l = left bearing 
r = rate constant 
I? = vertical 
x = x — direction 
y = y — direction 
bi = bias 

1 = top magnet 

2 = right magnet 

3 = bottom magnet 

4 = left magnet 

Superscripts 
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a = current amplifier 
c = controller 
l — lead 

p = sensor amplifier 
INTRODUCTION 

In recent years, the use of active magnetic bearings (AMB) for turbomachinery support 
has been an area of interest in both academia and in industry. Magnetic bearings 
provide the potential for significant improvements over other types of rotor supports, 
including elimination of wear, bearing friction-related energy losses and a means of actively 
suppressing rotor vibration. A critical feature of any magnetic bearing supported rotor 
system is auxiliary bearings to protect the soft iron core from rotor contact (and subsequent 
damage) during an overload or failure of the AMB. The present work is concerned with 
developing a suitable model and investigating the dynamic behavior of such a system. 

EXPERIMENTAL TEST FACILITY 

The model and control system development presented in this study is for an experi- 
mental test facility. The apparatus consists of a radial magnetic bearing that is supporting 
the right end of a rotor. The left end of the rotor is supported by a ball bearing suspended 
in a frame by four springs. A photograph and the corresponding schematic diagram are 
shown in Figure 1. 

The apparatus has two basic components: the rotor shaft and the magnetic bearing. 
The shaft is made of steel and is 0.0098 meters in diameter and 0.4572 meters in length. 
A steel disk of diameter 7.62 cm, thickness 2.54 cm, and mass 0.810 kg is placed at the 
midpoint of the bearing span. Threaded holes on the disk allow for imbalance to be added 
to the system. The rotor is driven by a variable speed motor with a controller. Shaft 
vibration is measured using eddy current proximity displacement sensors fixed to measure 
displacement in both the vertical and horizontal directions. 

The magnetic bearing is based upon a design described in Humphris, et al (1988). 
The basic parameters for the bearing are shown in Table 1. The bearing consists of 
four electromagnets equally spaced around a soft-iron core. A photograph and schematic 
diagram of the bearing are shown in Figure 2. 

SIMULATION MODEL 

A simulation model was developed for the rotor system described above. The model has 
three principal components - the rotor, magnetic bearing, and electronics. The governing 
equations for each are shown below. 

The rotor is modelled using the free-free bending mode shapes and natural frequencies 
obtained through finite element analysis. The finite element code uses 11 stations and the 
first four modes (two rigid body and two flexible modes) are included in the simulation 
model. The rotor equations of motion can be expressed in terms of modal coordinates as 


Qa: + firQy + W" Q x + ^ T F r — 0, 

(l.a) 

Qy-nrQ, + u;*Q 1 , + * T F y = 0, 

(1.6) 


The auxiliary bearing consists of a ball bearing with a clearance between the rotor and 
the stator. The governing equations are 
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where 


MqXq+C bxXq + KbxXb = Y x ,aux) ( l.c) 

Mb Yb+C By^B + Xb v Yb = F y y auxi (l.d) 

f x,lin(Nbl ) = ^ATn(^/a:Qa: "I" C/xQx) 

Fy l li n (Nbl') = *&Nbi ( A /t/Qy 4" C/j/Qj/) 

= ^l 2 ip6Cos6 + afasinO, 

F y ,i m b(N imb ) = FP^sinO - afocosO, 

F ar,aizx (-N& 3 ) = ^Ar 63 (FArcos/3 - F Rsin/3) 

F y y aux (-Y& 3 ) = $N b3 (FNsinf3 + Frcos/3) 

Fgrav = $ T MG 

F* = 4>K C (6 - D). 

If v r = 0, then F# = required static friction force to satisfy this condition. If v r ^ 0 
or if required static friction force > Fat, then = /^Fat. 

Fx = F x,lin H“ F x y imb “H F^a^a; Y Xyam b 


where 


• yjin + Fy^mb + Fj 


sin(/3) = 


Yh-Ys 

V'tXfi-Xs^ + tYfl-Ys) 2 , 


cos (/?) = 


Xfi-Xs 

v/(Xh-Xs ) 2 + (Yk-Y s ) 2 , 


with 


6 = yJ(*Q x -X b )* + (*Qy-Y b )*, 
cj) = 1 if 6 > A 
0 otherwise 

q x = <h T x r , 

Q y = * T Y r , 


X r ={X rl ,X r2 5***7 X rm } , 

Y r = {Y rl ,Y r 2 ,...,Y rm } T . 


The physical displacements of the rotor at the bearing and imbalance locations can be 
obtained using the following coordinate transformation: 
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(k — ^ 61 ^ 62 ^ 63 ? Nimb) 


N 

X r k = ^ ^ ki Qxii 

i—l 
N 

Yrk = ^ ^ ki Qyi ? 

i~l 

The position of the shaft at the magnetic bearing location is measured using horizontal 
and vertical proximity sensors. They are modelled as linear first order systems. 


d< = 

dt 

k v y - v p 

T P 

(2.a) 

dv p h 

k p x — vl 

(2.6) 

dt 

ft 

r p 


The active control system is based upon a proportional-derivative (PD) control law 
and implemented using analog circuitry. The controller is modelled as 


dv c v p 

-£ = k gV l>T' + (k a + k r )f t -T'v‘ v 


(3. a) 




(3.6) 


The control circuitry also includes a lead network. The governing equations are 


dVy 

dt 


= I7< + 


dyl 

dt 


7i v 'v 


(4.a) 


dv[ 

dt 


i _ 1 c . dvj 1 , 

t ri h+ dt tL h 


(4.6) 


Four current amplifiers (one for each coil) supply current to the magnetic bearing. 
They are modelled as 


djl 

dt 


~—ii + 


—I k{v[ 

j - a 1 v 


(5.a) 


di 2 
dt 


a„./ 


-71*2 + — ‘ k > 


(5.6) 
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(5.c) 


di 3 
dt 



dit 

dt 



(5 .d) 


The forces supplied by the magnetic bearing to the rotor are assumed to be decoupled 
in the horizontal and vertical directions. The effects of flux saturation are included in the 
model based upon the approach described by Lewis (1993). The equations for these forces 
are 


F x,amb — r 

( 2(D-x NB2 )+-f 

H-r 

i± ) 

( 2 (D + x N b2 ) + -^-) 2 

Hr 


(6.a) 


y,amb — /i<j.AlVj.( - 

( 2(D - y NB2 ) + 

Hr 

_____ i___) 

(2(D + y Ng2 ) + — ) 2 

Hr 


( 6 . 6 ) 


For design purposes, linearized models for the magnetic bearing forces are also 
needed. Linearization of (6.a) and (6.b) is performed about the nominal shaft position 
VNbi = Vn B 2 = 0. an< i x n b i = x n B 2 — 0 an d the nominal coil currents i\ = -i 3 = fcfvf* and 
*2 = — *4 = The effects of flux saturation are not included. The current stiffnesses 

that result are: 


K 


■i 

1 


H 0 Ak*v\' 'Ny 
2D> 


(7 -a) 


_ y-oAk^v^Nx 

" 2 ~ 2D 2 


(7.6) 


id = 


p 0 Aktv?N% 

2D 2 


(7.c) 
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Using the rotordynamic model described in the preceding paragraphs, a study of 
the combined dynamics of a magnetic bearing supported rotor with auxiliary bearings 
was performed. The parameters used in this work were identified from experimental 
evaluation of the test rig and are shown in Table 1. All of the responses illustrated in 
the figures are for this base parametric configuration, unless otherwise indicated. Figure 4 
shows the variation of system natural frequencies with rotor speed for the base parametric 
configuration (without auxiliary bearing contact). The rotor speed of 200 rad/sec was 
selected as a reasonable value near a critical speed in order to illustrate a worst case behavior 
of the system. 

Figure 5 shows the rotor response under normal operating conditions. The horizontal 
and vertical responses are of the same amplitude and frequency at steady state. The 
response consists of a transient region of about 0.5 seconds in length followed by steady 
synchronous oscillations. The initial conditions are set so that the rotor contacts the 
auxiliary bearing during the transient phase of the oscillations. The rotor vibrates 
downward, contacts the auxiliary bearing, bounces several times, and then settles into 
a steady synchronous whirl of about 0.05 mm in amplitude. 

Of particular interest in this study is the use of the remaining coils of the AMB (when 
one fails) to provide a sideloading force to reduce rotor vibration. The test configuration 
consists of allowing coil 2 to fail and investigating the effects of sideloading on the system. 
Coils 1 and 3 are kept in the same electronic configuration as if no failure had occurred. 
However, the voltage in coil 4 is set to a value that will pull the shaft horizontally over to 
the auxiliary bearing. Two basic scenarios for rotor control during failure are considered. 
Figures 6. a - 6.c illustrate the rotor behavior under the conditions for the first scenario, 
when the rotor speed is maintained at a constant value. For each case, the AMB fails at 
tj= 0.5 seconds. There is some perturbation of the vertical motion of the rotor when the 
horizontal coil fails due to coupling of the horizontal and vertical motion (primarily) through 
the gyroscopic terms. However, the steady state vertical responses are not significantly 
affected. Figures 7. a - 7.c show the rotor behavior for the same failure conditions (with 
sideloading) as above for the second scenario, when the rotor is decelerated at a constant rate 
upon AMB failure. Again the response amplitudes tend to be lowest for smaller clearance 
values. It is observed that the responses may be quite high as the rotor contacts the 
auxiliary bearing, even with sideloading from the remaining coil. However, sideloading is 
shown to be quite effective in encouraging rotor contact that will dissipate the vibrational 
energy. The response amplitudes tend to decrease as the auxiliary bearing clearance is 
decreased. However from a practical perspective, there are definite design trade-offs that 
must be considered. If the clearance is too small, the rotor may strike the auxiliary bearing 
during normal operation of the magnetic bearing if the rotor is perturbed. This tends to 
unnecessarily shorten the life of the auxiliary bearing. If the clearance is too large, the rotor 
vibration (upon AMB failure) will be excessively high and may damage the rotor or the 
magnetic bearing structure. 

The effect the sideloading voltage on the system dynamics is particularly interesting. 
The magnetic bearing forces are nonlinear functions of the dynamic clearance and the 
current. If one coil fails and the remaining coils are used to provide sideloading, the natural 
approach might seem to issue a constant command voltage, with the objective of having 
a constant sideloading force at steady state. However, the inherently unstable behavior of 
the magnetic bearing without active control coupled with the vibrational effects between 
the auxiliary bearing and the rotor produces a severely unstable response. A remedy is 
to command a voltage that is proportional to the horizontal shaft displacement. That is, 
v a = v“(D + x jv B2 ), where v* is a constant. With an ideal amplifier, any value of v“ can be 
used with such a system. However due to the delay in the dynamics from the time constant of 
the power amplifier, unstable behavior may also result if v% is too large. Figure 8 shows the 
maximum that can be commanded with stable behavior as a function of power amplifier 


8 



time constant. These results were developed by fitting a smooth curve through data points 
obtained by observing the behavior of the simulation model for various sideloading voltages 
and amplifier time constants. It must be concluded that power amplifier dynamics are a 
critical factor in the selection of appropriate levels of bearing sideloading. 

CONCLUSION 

A simulation model has been developed for a magnetic bearing supported flexible rotor 
system with auxiliary bearings including frictional effects. The model has been described 
in detail with experimentally obtained model parameters. The response characteristics 
for a variety of system configurations were presented and discussed, including the effects of 
sideloading. Some guidelines were given for the selection of appropriate levels of sideloading. 
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Figure l.b Schematic Diagram of Experimental Test Facility 






Figure 2. a Magnetic Bearing Assembly 



Figure 2.b Schematic Diagram of Magnetic Bearing 
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Figure 3 Block Diagram of Simulation Model 
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Figure 6. a Rotor Response With AMB Failure 
fi = 200rad/sec, Cl = 
tf — 0.5, oc — 0 
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Figure 6.b Rotor Response With AMB Failure 
Q = 200rad/sec, Cl = ^ 
tf = 0.5, a = 0 
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Figure 6.c Rotor Response With AMB Failure 
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Figure 7.a Rotor Response With AMB Failure 
Cl = 200rad/sec, Cl = 
if = 0.5, a = —100 
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Figure 7.b Rotor Response With AMB Failure 
f l = 200rad/sec, Cl, — 'rjr 
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Figure 7.c Rotor Response With AMB Failure 
Cl = 200rad/sec,Ci = y 
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Figure 8 Maximum Sideloading Force as 
a Function of Power Amplifier 
Time Constant, r a 
(for parameters of Table 1 ) 
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Abstract 

Foil bearings provide noncontacting rotor support 
through a number of thin metal strips attached around 
the circumference of a stator and separated from the 
rotor by a fluid film. The resulting support stiffness is 
dominated by the characteristics of the foils and is a 
nonlinear function of the rotor deflection. The present 
study is concerned with characterizing this nonlinear 
effect and investigating its influence on rotordynamical 
behavior. A finite element model is developed for an 
existing bearing, the force versus deflection relation 
characterized, and the dynamics of a sample rotor 
system are studied. Some conclusions are discussed 
with regard to appropriate ranges of operation for such 
a system. 

Nomenclature 

c xy = damping coefficient at x and y directions 
d = foil bearing clearance, mm 

f b = bearing restoring force acting on shaft, N 

I a = rotor polar mass inertia matrix, kg-m 1 2 

m = shaft and rigid disk mass, kg 

m e = imbalance mass, kg 

n f = number of foils 

N = total number of nodes on shaft 

N bI = node number at left bearing 

N b2 = node number at right bearing 

N disk = node number at disk on shaft 

p, 2 = foil bearing force coefficients 

q x = rotor modal coordinate vector in x 
direction, m 

q y = rotor modal coordinate vector 

in y direction, m 

r d = radius of rigid disk, mm 

r c = eccentricity of imbalance mass, mm 

Rj = journal radius, mm 

R b = bearing radius, mm 


t = foil thickness, mm 

x = displacement of shaft at x direction, m 

y = displacement of shaft at y direction, m 

Y = rotor free-free modal rotation matrix 

= rotor free-free modal displacement matrix 

r = Y T I a Y 

to = rotor speed, rad/s 

o> n = first critical speed, rad/s 

e = ratio of imbalance mass to total mass 

e bIb2 - bifurcation values of 8 

£ = damping factor 

cb n = matrix of rotor free-free 

natural frequencies, rad/s 
Mi,2 = mass modification coefficient 


Introduction 

The development of advanced types of high speed 
rotating machinery requires bearings that can withstand 
abusive environments for extended periods of time. 
Foil bearings offer a great deal of promise in this 
regard. Foil bearings have been applied successfully to 
wide range of light-duty, high-speed aerospace rotating 
machinery, auxiliary power units (APUs). They are 
presently employed in naval auxiliary power units 
(APUs) and numerous aircraft engines (e.g., Boeing 
747, 757, and 767, DC-10, F-15, F16, and Falcon 2000) 
[Heshmat (1994)]. 

Active magnetic bearing (AMB) supported rotor 
systems require auxiliary bearings to support the rotor 
and provide protection for the soft iron of the magnets 
in the event of AMB overload or failure. Bushings or 
rolling element bearings with a clearance between the 
support inner race and the rotor have typically been 
used for such applications. Recently, foil bearings have 
been considered for use as auxiliary bearings. They 
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offer the opportunity to provide continuous support to 
the rotor without significantly impacting the ability of 
the AMB to actively control rotor vibration (due to the 
inherent relative softness of foil bearings) while 
eliminating the need for oil lubrication and the 
associated cooling systems. Such a configuration could 
be quite advantageous for use in critical applications, 
such as jet engines, where safety is a primary 
consideration. 

There has been a great deal of research focussed on 
the dynamical behavior of foil bearings. Walowit, et al. 
(1973) first introduced a theoretical model to determine 
the static structural stiffness of foil bearings and 
performed an analysis of this model. In this analysis, 
the lubricant was assumed to be incompressible, the 
bearing was assumed to be an infinitely long, straight 
beam, and the friction forces due to foil/rotor contact 
were neglected. Bragin (1982) estimated the elastic 
characteristics of a foil sliding bearing with static 
loading using the foil packet model. Experimental 
results show that this idealized model makes it possible 
to determine the elastic reaction and stiffness of foil 
bearings. The foil packet, having low stiffness in 
comparison with the gas film, has the dominating effect 
on foil bearing elastic characteristics. 

Based upon the above discussion, an important 
area of concern is the expected dynamical behavior of 
a rotor supported by foil bearings. In order to 
investigate such a system, a flexible rotor simulation 
model is developed and analyzed using both numerical 
and analytical approaches. Some interesting results and 
observations are presented and discussed. 

Model Development 

The simulation model used in this study consists of 
foil bearings coupled with a flexible rotor. The model 
is described below. 

Foil Bearing Model 

The initial reference configuration chosen for the foil 
prior to the journal insertion in bearing is shown in 
Figure 1. This configuration is not unique and may 
differ in an assembled state from the free form state 
due to foil interaction loading. The relationships of the 
geometry parameters of foil bearing, for example, the 
number of foils n f , free foil length L f , foil thickness t, 
bearing radius R h , foil attachment angle <p, and foil 
radius of curvature R u are given in [Arakere and 
Nelson (1992 )]. The geometry of the foil bearing used 
in this study are listed in Table 1. 


Number of foils N f 

8 

Bearing Radius, R b 

45.95 mm 

Journal Radius, R^ 

40.00 mm 

Foil Radius, R, 

76.20 mm 

Bearing Length L 

106.7 mm 

Foil Thickness t 

0.229 mm 

Material 

Steel 


Table 1 : Geometric parameters and material 
properties of foil bearing 


Foil bearing clearance d 



Figure 1 Foil bearing configuration 

A finite element model of this foil bearing 
configuration is developed. The foils are modelled as 
elastic beams. There are several modelling issues that 
must be addressed in determining stiffness of foil 
bearing. These include how the foils interact, the impact 
of friction between the foils, and foil pre-loading. In 
order to address these problems, contact elements are 
defined between foils to simulate the mechanism of 
contact. Contact sliding and friction are considered in 
the analysis. To simulate the real case of shaft and 
bearings, the shaft is inserted into foil bearing, so that 
the foils are pre-loaded. The foil bearing force versus 
displacement characteristics, shown in Figure 2, are 
obtained by fitting a third order polynomial of the form 
(f p = PjX + p 2 x 2 3 ) to the finite element analysis output 
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data. This function is used to simplify the analysis. 
For the foil bearing used in the current study 
(parameters listed in Table 1), is 7970.0 N/m, and p 2 
is 1.496*10 lu N/m 3 . As can be seen from inspection of 
Figure 2, the approximation is quite good. 



Fiigure 2 Variation of foil bearing force with 
displacement 


flexible shaft 


rigid disk 

* 



Figure 3 Schematic diagram of simulation model 

The rotor shaft is modelled using the free-free 
bending mode shapes and natural frequencies obtained 
through finite element analysis. The finite element code 
uses 11 stations and the first four modes (two rigid 
body and two flexible modes) are included in the 
simulation model [Lawen and Flowers (1995)]. Based 
on these simplifications, the equations of motion for the 
system can be written as 


Flexible Rotor Model and Analysis 

The simulation study is based on the rotor/bearing 
configuration shown in Figure 3. It consists of a 
flexible rotor supported by foil bearings at both ends. 
A rigid disk with radius r d is placed at the midpoint of 
the bearing span. The rotor is supported by two foil 
bearings at each end. The "z" direction is the spin axis. 
The geometry and material property parameters used in 
this study are listed in Table 2. 


q + tor? + 5 y 2 q + <J> j F = 4> jT, (1) 

x y n x x x mb x 

q - v>Tq + 5 ) 2 q + <j / f = $ j f,. (2) 

y x n y y ymb v 

where 


Length of shaft 

620.0 mm 

Radius of shaft 

20.0 mm 

Radius of rigid disk 

100.0 mm 

Mass of disk 

2.0 kg 

Mass of shaft 

6.0 kg 

Material 

Steel 


Table 2 : Geometric parameters and material 
properties of rotor system 


(PiV p 2V c A )+ ®* 


(3) 


/=$ (P.9 +PA + c 4 )+® 

y 1 y 2 y y y 


N 

b2 
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the orbit of motion with large i/3 co subharmonic. In 

(6) both of cases shown in Figure 4-6 and Figure 7-9, the 
rotor system has exacdy the same parameters and rotor 
speed except the initial conditions. Comparison of 
vibration 

(7) 


( 8 ) 


(9) 

( 10 ) 



Figure 4 Purely synchronous time responses 


The physical displacements of the rotor at the right 
bearing, left bearing and rigid disk locations can be 
obtained by using coordinate transformations. 




A 


<") 


»„-EZ.*A <12) 

where M could be N hl , N h2 , or N disk . 

Numerical Simulation Analysis 

Numerical simulation is used to investigate the 
dynamical behavior of rotor system described above, 
e is taken as 0.4%, rotor speed co / co n is 4.93. 
Numerical simulation results, shown in Figure 4-9, 
show time response of the rigid disk in the horizontal 
and vertical directions, frequency domain analyses, and 
orbits. 

Figure 4-6 show that the vibration in time, frequency 
domain and orbit of motion. In this case, the system 
has purely synchronous vibration without subharmonic 
component. Figure 7 and 8 show that a large 1/3 co 
subharmonic response occurs besides the synchronous 
response in time and frequency domains. The 
amplitude of the subharmonic component, which is 
almost twice as large as the amplitude of the 
synchronous, dominates the vibration. Figure 9 shows 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
Frequency / rotor speed 

Figure 5 Purely synchronous frequency response 



X/d 

Figure 6 Purely synchronous trajectory 
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Figure 7 Time responses with a subharmonic 
component 



Figure 8 Frequency response with a subharmonic 
component 



Figure 9 Trajectory with a subharmonic 
component 


for these two cases shows that the rotor system 
supported by foil bearings has multiple response 
phenomena. It could be purely synchronous response 
or, it could have a combination of synchronous 
response and subharmonic component, depending upon 
the initial conditions. In the case of that system has 
purely synchronous vibration, the amplitude of vibration 
is much smaller than that of subharmonic case. 

Harmonic Balance Analysis 

The numerical harmonic balance approach is used to 
examine the synchronous and 1/3 <o subharmonic 
responses. The parameters are the same as in the 
previous study. Two case are examined : (1) the shaft 
is quite stiff (first two elastic natural frequencies are 
400 and 1200 hz) and (2) the shaft is quite flexible 
(first two elastic natural frequencies are 100 and 300 
hz). The results are shown in Figure 10. 



0 )/ 0 )„ 

Figure 10 Variation of response amplitude with 
rotor speed for harmonic balance analysis 

Figure 10 shows that the 1/3 co subharmonic 
components occur when rotor speed (o)/co n ) is larger 
then 3.0. The amplitude of subharmonic component is 
much larger than the amplitude of synchronous. It also 
shows that for the ranges of reasonable rotor speed 
c o/o) n (1.5— 8.5) and geometric configurations of shaft, 
the effect of rotor flexibility is not significant. It is 
concluded that the rigid rotor model described below 
adequately captures the dynamic behavior of such 
systems. In order to simplify the analysis and develop 
an understanding of the fundamental characteristics of 
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such systems, a rigid rotor model is now examined. 

Rigid Rotor Model and Analysis 

Based on the above discussions, the shaft is 
considered to be rigid. Furthermore, the masses of 
shaft and rigid disk are seem to be point masses. So, 
the rotor system has only two degrees of freedom of 
motion at vertical and horizontal directions. The bearing 
forces in the x and y direction are considered to be 
decoupled. So, the x and y equations are symmetric, 
and either of them can be selected to study without 
losing generality. The x equation is selected and is 
shown in (13). 


232 

jc + 2£<o x + co x + px = ec or cos(o)f) 


where 


(13) 


6pA 2 5*3cos(0 -20 )pAB 2 +3pB 3 +4Q 5=0 (18) 

12 3 


8 ZQB + 9 sin ( 0^20^45 = 0 


where 


2 2 2*2 
Q * co -o) , Q = co co , Q = a) — <o 
In 2 rt 3 « 9 


(19) 


( 20 ) 


Equations (16) - (19) have two possible sets of 
solutions - one for B=0 and the other for B*0. If B=0, 
equations (16) - (19) yield 


5=0 


( 21 ) 


co = 


2P. 


2P„ 


e = 


m 


(14) 


3 pA 3 + 4Q A = 4 er co 2 cos( 0 ) 


( 22 ) 


Harmonic Balance Analysis 

The harmonic balance approach is used to 
investigate the dynamical behaviors of the rigid rotor 
system for the range of operating speed (<o/o> n ~ 1. 5-8.5). 
Based on the results of the above simulations, it is 
assumed that the response, x(t), takes the form 


*(r) = Acos((*>t+QJ + Bcos(—oit+QJ 


(15) 


2 £Q A - er o> 2 sin ( 0 ) 

2 e 1 


(23) 


From (22) and (23), 0 can be canceled and an 
expression with only the synchronous response 
amplitude. A, the amplitude, can be obtained from 

(3 pA 3 + 4QA ) 2 + (2 5CM) 2 = 16e 2 r 2 co 4 (24) 


where A and B are amplitudes of the synchronous and 
l/3<o subhamonic amplitudes, and 0! and 0 2 are the 
phase shifts associated with the synchronous and 
subharmonic components. The effect of support 
damping is specifically included in this analysis. 
Substitution of (15) into (13) yields the following 
relationships between the synchronous and subharmonic 
response amplitudes and rotor speed. 

3pA 3 +6pAB 2 +cos(3 0^-0 ^ pB 3 +4 Q A =4 er co 2 cos(0 p 


If B*0, equations (25) - (28) can be simplified to yield 
3pA 3 +6pAB 2 +cos(3d -0 )pB 3 + 4Q A = 4 er <o 2 cos(0 ) 


2 r 


(25) 


8£Q 2 /1 + sin( 3Q^~Q^pB J = 4er co 2 sin ( 0p (26) 


(16) 6 ^ 2 + 3cos( 0 1 -29 2 )P/15 + 3 pB Z + 4Q 3 = 0 (27) 


85O 2 44 + sin(30 2 -0 i )pB 3 =4er co 2 siii(0 i ) 


(17) 


8£Q 2 + 9sin ( Q^ldJpAB - 0 


( 28 ) 
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Equations (24) - (28) show that a rotor system 
supported by foil bearings can have two possible 
responses. One is purely synchronous, and the other is 
a combination of synchronous and subsynchronous 
components. The analytical results demonstrate the 
correctness of conclusions obtained from numerical 
simulation of the flexible model. 

high synch, response 

unstable 

low synch, response 

synchronous comp. 

1/3 subharmonic comp. 



1 

x(t ) = /4cos(a)f) + 5cos(— ou) (29) 

Substitution of (29) into (13) yields 
3 pA 3 + 6pAB 2 + 3 pB 3 + 4Q A = 4er a> 2 (30) 

1 e v 7 


6pA 2 B + 3 pAB 2 + 3pB 3 + 40^2? = 0 (31) 

Here, A, B, p, Q, , and Q 3 have the same definitions 
as shown above. 

Bifurcation Analysis 

(1) Estimation of s ^ 

For estimation of s b2 , only the synchronous response 
exists. So, let B=0, and substitute B=0 into (30). It 
yields 

3pA 3 +4 Cl A = 4€r^co 2 (32) 


Figure 1 1 Variation of amplitude with imbalance 

Figure 1 1 shows that the rotor system supported by 
foil bearings could have multiple responses depending 
on the value of s . There are three possible responses, 
high synchronous, low synchronous, and a combination 
of synchronous and subsynchronous. If e< e bl , the 
system could have multiple synchronous responses 
(high or low), or a combination of 1/3 co 
subsynchronous and synchronous. The particular 
response the system actually takes depends on the 
initial conditions. If e bl < e< e b2 , the multiple response 
phenomena still exist (as in the situation for e< e bI ), 
except the subharmonic is of even order instead of odd 
order. Furthermore, if e> e b2 , the low synchronous 
response disappears. Determination of co b , e bl and 
e b2 are important for rotor designers. Estimation of 
those values could be done by using the necessary and 
sufficient condition for real solution(s) of the quadratic 
and cubic equation. Finally, the closed form of 
estimation of these values are developed in this paper. 

In order to explore these interesting phenomena, 
some assumptions must be made to simplify the 
analysis. First, damping is usually small, so it will be 
neglected. This results in 0,=O and 0 2 =O. Based upon 
these assumptions, (15) becomes 


Examination of (32) indicates that more than one real 
solution could exist if and only if A^O for (32). So, 


4er co 


4Q 


2 1 
3 ~ ) + 27 ( ^T 


) 3 < 0 


( 33 ) 


Assume, co > so Q=co 2 n - co 2 < 0. Solving (33), c b2 
is given 


e i — ( 1 
9 r 


2 P, 


m(l + |i)c» 


m( 1 + fi)<o - 2 p x 
2A 


( 34 ) 


where the total mass modification coefficient |J 2 due to 
the imbalance mass effect is 




2 

mu 


) 



( 35 ) 


(2) Estimation of e v _ 

Estimation of e M is very similar to the estimation 
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of e b2 , except that B*0. Under the present situation, the 
rotor system could be characterized by (30) and (31). 
Using B*0 condition, equation (31) can be simplified to 
yield 

6pA* + 3 pAB + 3 pB 2 + 4fl a = 0 (36) 


component could exist if and only if A > 0 of the B 
quadratic equation (36). So, co is given 


a) * 3 


2 

CO + 

71 


21A\ 

8 


- CO 


b 


(42) 


Because A and B « 1 in equation (30), the higher 
order terms of A and B can be neglected in order to 
obtain an estimation of A. So from equation (30), the 
estimation of A is given by if to * co n 

2 

er co 

« (T7\ 


Equation (42) shows that co b ^ 3co n , and a 1/3 co 
subharmonic could occur if co £ co b . This analysis result 
correlates well with the response observed in the 
flexible model numerical simulation. Comparisons of 
e bifurcation values (e M and e b2 ) to the corresponding 
values (obtained from numerical simulations of the 
flexible model) indicate that the relative error is in the 
range from 5% to 7%, which is quite reasonable. 


Substitution of equation (37) into equation (36) yields 
the quadratic equation in B 

2 2 
er o) er o) 

3 pB 2 * 3P-L—B + 6p( — — ) 2 + 4£1 3 = 0 (38) 

1 1 

In equation (38), B has real solutions if and only if A £ 
0. So, it is shown that 

2 2 
er o) er co 

A =9 p\ f - 16 p( 6p( — ^ ) 2 + 4Q 3 ),0 

1 1 


(39) 


Solving (39) for e, the bifurcation value e bl is given by 

2 P, 


4 2 

e* -( co 

3 r co 


(1 + u)m 


(1 + |r)mco 
42T 


(40) 


where the total mass modification coefficient p, due to 
imbalance mass effect is 


a 2 p 

4 , 2 1 V 

u -( co - — ) 

3r co m M 


mo) - 18/? i 
42p, 


(41) 


(3) Estimation of Ok 

Returning to equation (36), the subharmonic 


Discussions On Variations of p , and 

In the above analysis, the foil bearing parameters 
and p 2 ) are taken from the finite element analysis based 
the configuration listed in Table 2. It is of interest to 
consider the effect of variations of p, and p 2 on 
rotordynamical behavior. Figure 12 to Figure 15 show 
the higher, lower synchronous and 1/3 co subharmonic 
changes versus Pi and p 2 , where p 0l and p 02 are the 
values of p t and p 2 used in the previous analysis. 



Figure 12 Variation of amplitude with p 2 
(Purely synchronous) 



Figure 13 Variation of amplitude with p, 
(Purely synchronous) 
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P 2 /P(>2 “> 

Figure 14 Variation of amplitude with p 2 
(Subharmonic) 



Pi ! Poi 

Figure 15 Variation of amplitude with p ( 
(Subharmonic) 

Figures 12 and 13 indicate, that in the case of purely 
synchronous vibration, the higher synchronous 
responses decrease as p, and p 2 increase. The lower 
synchronous responses increase first and then disappear 
as pj and p 2 increase. Figure 14 and Figure 15 show 
the coexistence of synchronous and 1/3 co subharmonic, 
and indicate that the 1/3 co subharmonic components 
decrease and disappear as p t and p 2 increase. 
Furthermore, decreases of the subharmonic responses 
are much more sensitive to parametric variations than 
are the synchronous responses. 

Conclusions 

A study of the nonlinear dynamical behavior of a 
foil bearing supported flexible rotor system has been 
presented. The simulation and analysis work has been 


discussed in some detail. The following results and 
conclusions have been obtained. 

1. Very large 1/3 co subharmonic response are found 
if the rotor speed is approximately larger than three 
times of the first critical speed (co > <o b ). If the rotor 
speed is less than co h , the 1/3 co subharmonic response 
will disappear. Closed form expression for co b is 
developed. 

2. Multiple responses and bifurcation phenomena are 
shown. There are two particular values of imbalance 
mass ratio (e bl and e h2 ). If s < e bl , the system could 
have multiple synchronous responses (high amplitude 
and low amplitude) only if synchronous responses are 
excited. It could also have responses with a 
combination of low synchronous and 1/3 subharmonic 
responses. The particular response the system actually 
takes is determined by the initial conditions. If s bI < 
e < e b2 , the system may also have multiple 
synchronous responses (high amplitude and low 
amplitude) only if synchronous responses are excited, 
or it may have responses with a combination of low 
synchronous and even order subharmonic components 
instead of 1/3 o o components. If s > e b2 the system 
has only a single synchronous response. Also, closed 
form estimates of e bl and e b2 are developed. 

3. Variations of the linear and nonlinear coefficients of 
foil bearings have different impacts on the amplitudes 
of responses. If and p 2 increase, the high 
synchronous responses decrease. The low synchronous 
responses are not affected significantly. In the case of 
the existence of 1/3 co subharmonic responses, the 
subsynchronous responses decrease much faster than 
high synchronous responses. 
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ABSTRACT 


This study investigates the application of syn- 
hronous interaction dynamics methodology to the 
lesign of auxiliary bearing systems. The technique is 
applied to a flexible rotor system and comparisons are 
nade between the behavior predicted by this analysis 
nethod and the observed simulation response char- 
acteristics. Of particular interest is the influence of 
oupled shaft /bearing vibration modes on rotordy- 
Lamical behavior. Experimental studies are also per- 
ormed to validate the simulation results and provide 
asight into the expected behavior of such a system. 


NOMENCLATURE 

a = rotor polar mass inertia matrix, kg-m 2 

( = stiffness, N/m 

A b = auxiliary bearing mass, kg. 

i = total number of modes considered 

f b i = node number at leftmost bearing 

^62 = node number at rightmost bearing 

= node number at auxiliary clearance bearing 
} x = rotor modal coordinate vector in X direction 
) y = rotor modal coordinate vector in Y direction 
= time, s 

l b = amplitude of bearing vibration, m 
l r = amplitude of rotor vibration, m 
f b = auxiliary bearing physical coordinate vector 
in X direction, m 


* Work performed while at Auburn University 


Y b = auxiliary bearing physical coordinate vector 
in Y direction, m 

X r = rotor physical coordinate vector in X 
direction, m 

Y r = rotor physical coordinate vector in Y 
direction, m 

V = structural damping coefficient, N/m 
W = linear damping coefficient, N-s/m 
A = radial clearance in auxiliary bearing, m 
= rotor free-free modal rotation matrix 
$ r = rotor free-free modal displacement matrix 

r = 

f l = rotor operating speed, rad/s 
u> n = matrix of rotor free-free natural 
frequencies, rad/s 
if) = imbalance vector 

Subscripts 

bx = auxiliary bearing, x-direction 
by = auxiliary bearing, y-direction 
c = contact 

xr = rightmost bearing, x-direction 
xl = leftmost bearing, x-direction 
yr = rightmost bearing, y-direction 
yl — leftmost bearing, y-direction 

INTRODUCTION 

Auxiliary bearings for magnetic bearing sup- 
ported rotor systems serve two basic missions. First, 
they protect the soft iron of the magnetic bearings 
during speed runup (through critical speeds) and 
overload. Second, they provide added rotor support 
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luring critical operating conditions. A vital require- 
nent is that the auxiliary bearings should perform 
heir roles without interfering with the rotor during 
lormal operation of the system. 

There has been quite a bit of work over the past 
lecades concerned with the dynamics of a rotor in- 
fracting with a clearance bearing. Some of the most 
>rominent work in the area has been performed by 
£hrich (1965, 1966, 1988), Bently (1974), Muszyn- 
ka (1984), Childs (1979 and 1982), Choi and Noah 
1987), and Kim and Noah (1989). The work of 
Cirk and Ishii (1993) figures prominently in the liter- 
Lture on this subject, with regard to auxiliary bear- 
ng applications. Schweitzer, et al. (1994) presents 
t, good review and discussion of issues related to the 
ouch-down dynamics of rotors on auxiliary bearings. 
,awen and Flowers (1995) provides a more detailed 
eview of the relevant past literature on the subject. 

While these studies have greatly enhanced the 
mderstanding of the dynamics of rotors supported 
\y auxiliary clearance bearings, there is still a need 
o develop a better understanding of the possible 
[ynamics of such systems and to devise guidelines 
or the design of auxiliary bearings. With this idea 
a mind, motivation for the current work is due to 
, study conducted by Black (1968). He examined 
he synchronous interaction of a rotor and stator 
ue to contact across a clearance annulus. By 
efining the rotor and stator displacement s,f r and 
5 respectively, in terms of their polar receptances 
or transfer functions), the displacements can be 
xpressed in terms of the forces acting in the system 
s shown below: 


= —an Pc + \pW 9 

(l.a) 

r a = flu P c 

(1.6) 


where an and /?n are the polar receptances of 
be rotor and stator respectively, P c is the force due to 
ontact, and | p\e j6 is the whirl radius of the rotor in 
be absence of contact due to imbalance. The radial 
learance,12 , can be defined in terms of the stator 
isplacement relative to the rotor during interaction 
s shown in (2). 


R = r r -r s (2) 

Substituting (l.a) and (l.b) into (2) yields 


R + \a n +flu\Pe jl) = \p\e 36 (3) 


where T) is the phase lead of (an -f- flu) relative 
to P as shown in (4). 

(an + flu) = | <*n + Pu\e 3V (4) 

Obviously if the rotor whirl radius due to imbal- 
ance exceeds the radial clearance, \p\ > |fZ| , then in- 
teraction occurs between the rotor and the stator due 
to contact. However, using vector triangles, Black 
also finds a criteria where interaction is possible even 
if \p\ < |i2| as shown in (5). 


P -cos(77)-y^|y-sin 2 (»7) 

\p\ = |an+0u| (5) 

For contact to occur, P, the contact force, must 
be positive which yields the following two conditions: 


t] < —90 (6.a) 


I s Hv)\ < (6.6) 

These equations can be used to predict when 
interaction between stator and rotor is possible for 

\ P \ < \R\. 

In essence, by defining the rotor and stator 
(bearing or bearing/housing) in polar receptance form 
and determining the phase lead, 7/, associated with 
the combined polar receptance, zones of interaction 
can be defined where rotor and stator coupling can 
occur even for rotor imbalance responses smaller than 
the radial clearance. 

The present paper applies the synchronous inter- 
action methodology devised by Black (1968) to the 
problem of auxiliary bearing design. The method is 
applied to a rotor and housing system modelled us- 
ing modal coordinates for the structure, simulation 
and experimental studies are described, and the re- 
sults compared and discussed. The implications for 
auxiliary bearing design are detailed. 

ANALYSIS 

Figure 1 illustrates the basic simulation model 
used for this investigation. It consists of a flexible 
rotor supported at both ends by magnetic bearings. 
A rigid disk with adjustable imbalance is placed at the 
midpoint of the bearing span. As a simplification, the 
magnetic bearings are modeled as spring and dashpot 
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:ystems. Interaction with only a single auxiliary 
>earing is considered. 

The rotor is modeled using the free-free bend- 
ng mode shapes and natural frequencies obtained 
hrough finite element analysis. Using this modal 
lata, the equations of motion for the flexible rotor 
nodel can be developed, as described in Ryan (1991). 

Figure 2 shows the auxiliary bearing model used 

0 study the effects of the bearing properties on 
he rotor response during coupled behavior between 
he flexible rotor and the auxiliary bearing. As 

1 simplification, the auxiliary bearing was modeled 
is an antifriction bearing with a clearance. The 
auxiliary bearing housing is represented as a mass 
attached with springs to the bearing. 

Using the methodology described above, poten- 
ial zones of coupling between the flexible rotor and 
he auxiliary bearing and between the flexible rotor, 
he auxiliary bearing, and the auxiliary bearing hous- 
ag are determined. The basic approach is to deter- 
nine the polar receptance of the rotor and the stator 
auxiliary bearing or auxiliary bearing/housing) and 
hen find the phase lead of the combined rotor-stator 
tolar receptance. If the phase leads meets the crite- 
ia of (6.a) and (6.b), then interaction between the 
iexible rotor and stator is possible. 

With this in mind, the equations of motion 
,re written in complex notation using the following 
oordinates: 


Qr — Qx 4 ' jQy 

(7.a) 

Rb = Xb + jYb 

(7.6) 


The stiffness and damping terms for a particular 
ystem component, e.g. the left magnetic bearing, 
/ere also equated in the horizontal and vertical 
irections, e.g. Kt = K x i = K y i. The resulting 
quations of motion are: 

Q r -jttTQ r + U 2 r Qr 
+*T*rjni(KlQr + (f + W ‘)0r) 

+* T *r, Nb 2(K r Q r + (~ + W r )Q r ) 

+$jF nl = $jF rtimb (8 .a) 

Mb&b + (^ + W b )kb + K b R b = F nl (8.6) 


F n i and F, m j, refer to the contacting and imbal- 
ance forces respectively. 

Assuming a synchronous solution for the rotor 
and bearing, 

Qr = Q ro e jClt (lO.a) 

R b = R bo e jUt , (10.6) 

substitution of (10. a) and (lO.b) into (8.a) and 
(8.b), yields the following polar receptances for the 
rotor-auxiliary bearing model: 

a n = [^n — + T) + $rF S upp]~~ X (ll.a) 

fill = [Kb - M b tf + jft(^ + Wb)]- 1 (11.6) 

where 


F aupp = ,mi[K, + + W,)] 

+ $r,/V62[tfr + + W r ) ] (ll.c) 

The receptances, an and /3 n , reflect the equa- 
tions of motion for the rotor and the auxiliary bearing 
being written in the form below. 

Qr = a n $J(F imb - F nt ) (12.a) 

Rb = Pll+Fni (12.6) 

It should be noted that the polar receptance 
given by (11. a) is in terms of the rotor modal co- 
ordinates. For the receptance in rotor physical co- 
ordinates at the auxiliary bearing location, 3, the 
following transformation must be performed. 


N 

3 a ll (13) 

*=1 

The combined rotor-bearing polar receptance can 
then be written as follows: 

ar + 0n = \a r +0 n \e jr) (14) 

where 77 is the phase lead of interest. 

Therefore, (ll.a), (ll.b), and (13) are combined 
to form the combined polar receptance (14) for the 
auxiliary bearing interacting with the flexible rotor. 
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Rigid Disk 


Auxiliary Bearing 



Figure 1 Schematic Diagram of Simulation Model 

Rotor 


auxiliary clearance bearing rig. The shaft is made of 
steel and is 0.374 inches in diameter and 18.0 inches 
in length. It is supported at 1.0 inch from the right 
end by bail bearings suspended in a frame bv four 
springs and at 1.0 inch from the left end by a bushing 
with a tight clearance. These supports represent the 
magnetic bearings. The stiffness of the left support 
is 17511 N/m. for both the horizontal and vertical 
directions. It is used to somewhat isolate the rotor 
from the effects of the flexible coupling which attaches 
the rotor to the motor and to enforce low amplitude 
vibration at this location to protect the motor. The 
stiffness of the right support, for both horizontal and 
vertical directions, is 2539 N/m. This lower stiffness 
allows for significant vibration of the rotor in the 
speed range of the motor (0-10,000 rpm). A rigid 
disk with holes for placing imbalance screws is placed 
at the midpoint of the bearing span. 



The phase lead, 77 , associated with the combined 
:omplex polar receptances is then found over a range 
}f operating speeds. If the conditions of (6. a) and 
6.b) are met for rotor imbalance responses less than 
uhe clearance, then interaction between the rotor and 
auxiliary bearing is possible. 

EXPERIMENTAL MODEL 

Experimental tests were performed in order to 
validate the behavior predicted by the simulation 
model and to gain some insight into the dynamical 
responses that are to be expected. A photograph of 
the experimental test rig is shown in Figure 3. It 
has two basic components: a flexible shaft and an 



Figure 3 Photograph of Experimental Apparatus 


The rotor is driven by an adjustable speed mo- 
tor with feedback speed controller. Shaft vibration is 
measured using eddy current proximity displacement 
sensors fixed to measure displacement in the vertical 
and horizontal directions. During testing, the dis- 
placement signals of the shaft are sent to the signal 
analyzer where the frequency and amplitude compo- 
nents of the shaft response are recorded. The results 
are compared with those from the analytical stud- 
ies. The auxiliary clearance bearing/housing assem- 
bly consists of a bushing suspended in a frame by four 
springs and is situated at the right end of the rotor. 

DISCUSSION 

In order to evaluate the effectiveness of the anal- 
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rsis technique described above, a series of parame- 
er variation studies were performed. The nominal 
jarametric configuration is shown in Table 1 and all 
he cases presented use these parameters unless other- 
vise indicated. These parameters are representative 
>f reasonable values for the experimental test facility. 
The first four rotor modes (two rigid body and two 
lexible modes) are included in the simulation model. 

The dynamic behavior of such a system is gov- 
erned by shaft flexibility, auxiliary bearing stiffness, 
lamping, and mass, imbalance, auxiliary bearing lo- 
ations, bearing stiffness and damping, and auxiliary 
rearing clearance. The parameters over which the 
rearing designer has the most control are the auxil- 
ary bearing stiffness and the clearance, which are the 
ocus of the discussion below. 

First, let us begin by considering the possible 
>ehaviors for different auxiliary bearing clearances, 
figures 4.a - 4.e show some representative results for 
everal values of bearing clearance. Figure 4.a shows 
he variation of 77 with the rotor speed. (Naturally, it 
3 independent of the values of bearing clearance that 
,re chosen.) Similarly, Figure 4.b shows the varia- 

ion of ^ and sin 77 with rotor speed for the values of 

uxiliary bearing stiffness that are considered. Com- 
parison of Figure 4.a with Figure 4.b indicate that the 
peaks in the 77 curve tend to occur near the respective 
ombined rotor/stator critical speeds, as indicated by 

he amplitude peaks of the y and sin 77 curve in Fig- 
ure 4.b. For rotor speeds oelow the first critical, 
ateraction is not possible. For rotor speeds above 
he third critical, 77 approaches zero and interaction 
3 possible for all rotor speeds for which sin 77 > 

/hich is the entire remaining range of rotor speed 
hown. For each of the three values of clearance, there 
3 a region between the first and second critical for 
fhich steady rotor/bearing interaction is impossible 
l the rotor imbalance response amplitude, p, is lower 
han the clearance. As the clearance is increased, 
his regime expands dramatically from a quite narrow 
and (196 to 247 rad/sec) to a rather broad range of 
otor speeds (0 to 430 rad/sec). Inspection of Fig- 
res 4.c (rotor response amplitudes) and 4.d (bearing 
esponse amplitudes) reveals some interesting trends, 
'igure 4.d shows that interaction can occur over vir- 
ually the entire range of rotor speeds (except for the 
egime between the first and second critical speed as 
escribed above), which validates the behavior pre- 
icted using Black’s interaction criteria. At the three 
ritical speeds, the auxiliary bearings tend to attenu- 
te the peak response amplitudes for all cases. Very 
ttle shift is seen in the location of the first criti- 
al speed as a result of interaction with the auxiliary 
earing. These results agree well with experimental 
at a for the region 0 < 12 < 1000 rad/sec, as shown 


in Figure 4.e. 

Next, let us examine the effects of variations in 
auxiliary bearing stiffness on possible interactions be- 
tween the rotor and the bearing. Some results for a 
variety of auxiliary bearing stiffnesses are summarized 
in Figures 5. a - 5.e. Figure 5. a shows the variation in 

77 and Figure 5.b shows and sin 77 for the range of 

rotor speed and the values of auxiliary bearing stiff- 
ness that are considered. Again, no interaction is pos- 
sible at very low rotor speeds (12 < 50) and interaction 
can occur for all rotor speeds above the third criti- 
cal speed (fi > 1225). In between these extremes, 
a variety of behaviors are possible depending upon 
the rotor speeds. In contrast with Figure 4 .b, Figure 

5.b indicates that sin 77 < ^ for almost the entire 


range of rotor speed. Unlike the previous case (for 
varying clearance values), the interaction dynamics 
here are dominated by the requirement that 77 < —90 
degrees. Inspection of Figure 5.a indicates some in- 
teresting trends. For the lowest value of auxiliary 
bearing stiffness (Jf &= 88 N/m), interaction is possible 
over the entire range of rotor speed. As the auxiliary 
bearing stiffness is increased, there is again a range 
of rotor speed between the first and second critical 
speeds for which interaction is not possible. For the 
highest value of auxiliary bearing stiffness (^=2625 
N/m), interaction is not possible for rotor speeds be- 
tween 144 and 401 rad/sec. Figures 5.c (rotor re- 
sponse amplitudes) and 5.d (bearing response ampli- 
tudes) again validate the analysis results predicted 
using Black’s synchronous interaction criteria. The 
auxiliary bearings attenuate the peak response ampli- 
tudes for all critical speeds. The first combined criti- 
cal speed is increased by the presence of the auxiliary 
bearings (indicating dominance of the stiffness contri- 
bution from the auxiliary bearing), while the second 
and third combined critical speeds are reduced (indi- 
cating dominance of the mass contribution from the 
auxiliary bearing). These results are further validated 
by the experimental data for the region 0 < SI < 1000 
rad/sec, as shown in Figure 5.e. 


CONCLUSIONS 

Typically bearing support structures are de- 
signed considering only the stiffness and damping 
characteristics of the structure. For standard appli- 
cations, this is an acceptable approach. However, the 
usage of auxiliary clearance bearings introduces an 
additional parameter, clearance, and an essentially 
nonlinear feature into the rot or dynamics. In order 
to properly design such systems, it is necessary to 
develop a detailed understanding of the sorts of dy- 
namical behavior that may be possible. Toward this 
end, the method of Black for the analysis of inter- 
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2.0 


Parameter 

Value 

Units 

Rbxi K-by 

471 

N/m 

K c 

87,557 

N/m 

K lx , K, y 

17,510 

N/m 

Krxi K f y 

2,539 

N/m 

M b 

0.033 

kg 

Vbx, V by 

2,000 

N/m 

v,x, V, y 

1,000 

N/m 

V rx , v ry 

600 

N/m 

Wbx, W by 

0 

N-s/m 

w.x, w ly 

0 

N-s/m 

W rx , w ry 

0 

N-s/m 

4>i 

3.0 x 10" 6 

kg-m 

i’w 

2.5 x 10- 4 

kg-m 

A 

5.0xl0 -5 

m 


Table 1 Simulation Model Parameters 



Figure 4.a Phase Lead for Varying 

Auxiliary Bearing Clearances 

f 

kting rotor/stator systems has been applied to the 
fase of a flexible rotor-bearing system, where the ro- 
jor is modelled using modal coordinates. A simple 
lexible rotor system is studied, with simulation and 
Experimental studies being conducted for a variety of 
parametric configurations. These results have been 
presented and discussed. Regions of interaction have 
been detailed and the overall effects of various auxil- 
ary bearing clearance and stiffness values have been 
lescribed. Provided that the assumption of purely 



Rotor Speed (rad/sec) 


Figure 4.b Noncontacting Imbalance Response 
Amplitudes and sin(n) for Varying 
Auxiliary Bearing Clearances 



Rotor Speed (rad/sec) 

Figure 4.c Simulation Rotor Response 
Amplitudes for Varying 
Auxiliary Bearing Clearances 

synchronous vibration is reasonably satisfied, Black’s 
method for identifying regions of synchronous inter- 
action is shown to be a versatile tool for the analysis 
of rotor systems with auxiliary clearance bearings. 
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Rotor Speed (rad/sec) 


Figure 4.d Simulation Bearing Response 
Amplitudes for Varying 
Auxiliary Bearing Clearances 



Rotor Speed (rad/sec) 

Figure 4.e Experimental Rotor Response 
Amplitudes for Varying 
Auxiliary Bearing Clearances 
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Figure 5.b Noncontacting Imbalance Response 
Amplitudes and sin (77) for Varying 
Auxiliary Bearing Stiffness 

REFERENCES 

Black, H. F., 1968, “Interaction of a Whirling 
Rotor With a Vibrating Stator Across a Clear- 
ance Annulus,” Journal Mechanical Engineering 
Science, Vol. 10, No 1, pp. 1-12. 

Childs, D. W., 1982, “Fractional-Frequency Ro- 
tor Motion Due to Nonsymmetric Clearance Ef- 


7 







R b (meters) R r (meters) 



0 500 1000 1500 2000 


Rotor Speed (rad/sec) 

Figure 5.c Simulation Rotor Response 
Amplitudes for Varying 
Auxiliary Bearing Stiffness 



0 500 1000 1500 2000 

Rotor Speed (rad/sec) 

Figure 5.d Simulation Bearing Response 
Amplitudes for Varying 
Auxiliary Bearing Stiffness 

fects,” ASME Journal of Engineering for Power, 
Vol. 104, pp. 533-541. 

Ehrich, F. F., 1965, “Bistable Vibrations of 
Rotors in Bearing Clearance,” ASME Paper 65- 
WA/MD-1. 

Ehrich, F. F., 1988, “High Order Subharmonic 



0 200 400 600 800 1000 


Rotor Speed (rad/sec) 

Figure 5.e Experimental Rotor Response 
Amplitudes for Varying 
Auxiliary Bearing Stiffness 

Response of High Speed Rotors in Bearing Clear- 
ance,” ASME Journal of Vibration, Acoustics, 
Stress, and Reliability in Design, Vol. 110, pp. 
9-16. 

Ishii, T., and Kirk, R. G., 1991, “Transient Re- 
sponse Technique Applied to Active Magnetic 
Bearing Machinery During Rotor Drop,” DE- 
Vol. 35, Rotating Machinery and Vehicle Dy- 
namics, ASME, pp. 191-199. 

Kim, Y. B., and Noah, S. T., 1991, “Steady-State 
Analysis of a Nonlinear Rotor-Housing System,” 
Journal of Engineering for Gas Turbines and 
Power, Vol. 113, pp. 550-556. 

Lawen, James L., Jr., and Flowers, George 
T., 1995, “Synchronous Dynamics of a Coupled 
Shaft/ Bearing/Housing System With Auxiliary 
Support From a Clearance Bearing: Analysis 
and Experiment,” presented at the ASME Inter- 
national Gas Turbine Congress and Exposition, 
Houston, Texas, June 5-8, 1995, ASME paper 
no. 95-GT-216. 


8 






